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ABSTRACT 


This  thesis  deals  with  the  pore  pressure  behaviour  within 
embankment  dams  resulting  from  resevoir  impounding  and  the 
transition  to  steady  state  seepage.  The  theories*  analyses 
and  results  are  specifically  related  to  the  case  history  of 
the  Mica  dam  in  British  Columbia. 

The  processes  which  affect  the  pore  pressures  within 
the  core  during  resevoir  loading  and  transient  seepage  have 
been  discused  and  analysed.  The  relative  contribution  of 
each  mechanism;  total  stress  change*  transient  seepage  and 
consolidation,  has  been  evaluated.  As  a  product  of  this 
analysis  hydraulic  fracture  has  been  investigated  and  some 
conclusions  are  drawn  regarding  its  likelihood  •  Steady 
state  seepage  analyses  have  been  carried  out  for  the  Mica 
dam,  and  the  influence  of  such  parameters  as  anisotropy* 
non-homogeneity  and  underseepage  have  been  investigated.  A 
short  discussion  on  piezometer  reliability  is  included 
within  the  text. 

The  effect  of  ins trument— lead  placement  technique  has 
been  investigated  both  in  a  general  sense  and  specifically 
related  to  the  Mica  dam. 

It  is  concluded  as  a  result  of  these  studies  that  the 
method  of  instrument  installation  is  critical  for  narrow 
core  dams.  The  method  employed  in  the  Mica  dam  has  caused 
local  anomalies  to  the  pore  pressure  distribution  and  so  the 
core  piezometers  are  not  reflecting  the  average  pore 
pressure  distribution. 
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CHAPTER  1 


INTRODUCTION 


1-a.l  Mgdegp  Jag  EggjgQ 

During  the  last  two  decades  there  have  been  significant 
advances  in  embankment  design  techniques*  An  improved 
understanding  of  the  physical  processes  occurring  during 
embankment  construction  has  been  achieved  largely  through 
the  development  of  analytical  methods  for  stress  and  strain 
analysis  coupled  with  an  increase  in  the  use  of 
instrumentation  to  monitor  embankment  performance. 

Although  present  day  analytical  tools  are  powerful  and 
sophisticated,  the  designers  predictive  capabilities  are 
reduced  by  an  incomplete  knowledge  of  the  constitutive 
relationships  for,  and  interaction  between,  the  composite 
soil  types.  These  laws  and  interactive  behaviour  mechanisms 
are  being  refined  through  careful  b ack— ana lys i s  of  well 
instrumented  case  histories. 

In  current  practice  the  finite  element  method  (F.E.M. ) 
is  often  used  to  provide  insight  into  possible  problem  areas 
in  the  design.  It  may  also  provide  a  rough  qualitative 
assesment  of  the  expected  Instrument  readings.  The  F.E.M.  is 


1 
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used  for  the  analysis  of  stress*  deformation*  consolidation* 
and  seepage  conditions  within  the  embankment* 

The  factor  of  safety  against  slope  failure  is 
determined  by  the  use  of  conventional  limit  equilibrium 
methods*  These  methods  employ  the  principle  of  effective 
stress  which  requires  the  prediction  of  the  pore  pressure 
distribution  within  the  embankment*  Conventionally  factors 
of  safety  are  calculated  for  the  end  of  construction*  first 
filling,  steady  state  and  rapid  drawdown  conditions  in  order 
to  establish  the  critical  case  when  the  fcctor  of  safety  is 
a  minimum*  The  pore  pressure  distribution  for  both  the  end 
of  construction  and  steady  state  cases  are  well  defined  by 
routine  analysis*  however  the  first  filling  and  rapid 
drawdown  conditions  require  simplifying  assumptions  for  the 
pore  pressure  distributions  in  order  to  continue  with  the 
ana lys is* 

It  is  apparent  that  little  attention  has  been  given  to 
estimating  the  pore  pressure  distribution  during  the  period 
between  the  end  of  construction  and  the  establishment  of 
steady  state  seepage*  which  is  defined  here  as  the 
'transient*  seepage  stage*  The  reason  for  this  is  the 
complex  interaction  of  the  processes  involved  during 
resevoir  loading*  Clearly  it  is  desirable  to  gain  further 
insight  into  the  behaviour  of  the  pore  pressure  distribution 
within  the  embankment  during  and  subsequent  to  resevoir 
filling*  This  is  particularly  important  because  the  first 
filling  of  a  resevoir  is  considered  as  a  critical  stage  in 
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the  working  history  of  a  dam* 

A  recent  addition  to  modern  dam  design  is  the 
evaluation  of  hydraulic  fracture  potential*  This  is  a  poorly 
understood  phenomenon  which  certainly  requires  further 
research  and  clarification* 

Even  though  the  embankment  dam  designers*  technical 
expertise  and  understanding  is  improving  rapidly  there 
remain  periods  in  the  life  of  an  embankment  which  are  not 
well  understood*  The  most  critical  of  these  is  the 
'transient'  seepage  stage  •  There  may  be  significant  changes 
in  the  pore  pressure  distribution  during  this  phase  which 
the  designer  is  at  present  unable  to  analyze  accurataly*  The 
early  part  of  this  period,  the  resevoir  filling  stage,  is 
probably  the  most  critical  period  for  the  eveluation  of 
hydraulic  fracture  potential,  which  may  be  of  major  concern 
with  some  embankment  configurations* 

Aim  slI  Ihls.  Ihesis 

This  thesis  was  undertaken  with  the  following 
objectives : 

1«  To  take  a  carefully  instrumented  and  well 

documented  case  history,  namely  the  Mica  Dam, 
and  investigate  the  pore  pressure  distribution 
within  the  embankment  during  the  'transient' 
seepage  stage  (between  end-oi-cons  true  ti on  and 
steady  seepage)* 

2*  To  define  the  processes  taking  place  during  this 


, 


* 


4 


in 

du 

re 


period  which  cause  alterations  in 
pressure  distribution,  and  to  eval 
relative  contribution  of  each  proc 
3*  To  attempt  to  evaluate  the  hydraul 
potential  in  the  core  of  this  dam 
resevoir  filling* 

4*  To  briefly  discuss  the  performance 
pressure  measuring  devices  install 
In  summary  this  thesis  is  an  attempt 
sight  into  the  behaviour  of  the  pore  pre 
ring  and  subsequent  to  resevoir  filling, 
ference  to  the  Mica  Dam* 


the  pore 
uate  the 
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ic  fracture 
during 
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to  provide 
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with 

spec 
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1_»  2  Qalllns  Q.1  ths  i.frggls 

The  thesis  is  organised  in  the  following  sequence: 
Chapter  2  provides  a  general  description  of  the  Mica 
Dam  including  a  discussion  on  the  piezometer  types  and 
methods  of  installation* 

Chapter  3  discusses  the  reliability  and  performance  of 
the  piezometers  in  general* 

Chapter  4  details  the  steady  state  seepage  analysis* 
Chapter  5  compares  the  observed  and  calculated  pore 
presure  distribution  within  the  core* 

Chapter  6  analyses  the  seperate  physical  processes 
involved  during  resevoir  filling*  It  studies  the  transient 
flow  and  evaluates  the  hydraulic  fracture  potential  within 


the  core*  Finally  these  results  are  compared  to  the  observed 
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pore  pressure  distribution* 

Chapter  7  evaluates  the  effect  of  imperfect  instrument 
installation  procedure  both  generally!  and  specifically  for 
the  Mica  Dam  core*  These  results  are  compared  to  the 
observed  distribution  calculated  from  the  piezometer 
readings. 

Chapter  8  presents  the  conclusions  of  the  thesis  and 
discusses  some  recommendations  for  further  research* 


CHAPTER  2 


THE  MICA  DAM 


2.i  1  general  I2ggg.rjtp.tlgja 

Mica  dam  is  situated  on  the  Columbia  River  in  British 
Columbia  immediately  west  of  the  Rocky  Mountains*  and  135 
kilometers  north  of  Revelstoke*  The  site  is  in  a  narrow 
gorge  where  the  Columbia  River  leaves  the  Rocky  Mountain 
Trench#  The  location  is  shown  in  Fig*2*l* 

The  dam  was  constructed  as  part  of  the  Columbia  River 
Treaty  (1964)*  which  provides  regulation  for  maximum 
production  of  power  and  flood  control  of  the  river*  The 
project  includes  a  2500  MW  powerplant  which  develops  the 
head  created  by  the  dam*  The  storage  provided  by  the  dam  is 
approximately  25  x  109  m3  and  the  length  of  the  reservoir 
created  is  about  208  kilometers*  from  Golden  to  Valemont  in 
the  Trenc h • 

The  dam  is  a  zoned  earthfill  embankment*  with  a 
near— vertical  central  impervious  core  of  heavily  compacted 
glacial  till*  The  shells  are  compacted  sand  and  gravel  of 
which  the  outer  edges  are  compacted  and  dumped  rock*  The 
general  fill  arrangement  is  shown  in  Figure  2*2*  The 
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Location  of  Mica  Dam 
(after  Law  1975) 


FIGURE  2.1 
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(1)  Axis  of  dam 

(2)  Upstream  cofferdam. 

(3)  Downstream  cofferdam. 

(4)  Core-glacial  till. 

(J)  Shells-sand  and  gravel. 

(6)  Outer  shells-sand  and  gravel  or  rock. 

(7)  Drawdown  zone-gravel,  cobbles,  boul¬ 
ders  or  rock. 

(8)  Slope  protection. 

(9)  River  overburden. 

(10)  Bedrock. 

(11)  Original  riverbed. 

(12)  Limit  of  blanket  grouting. 

(13)  Deep  exploratory  holes. 

(14)  Normal  maximum  reservoir  level 

(15)  Minimum  drawdown  level 


FIGURE  22  Mica  Dam  -  General  Fill  Arrangement 

(after  Nussbaum  1976) 


embankment  was  constructed  over  4  seasons  between  April  1969 


and  November  1972*  The  reservoir  filling  commenced  in  April 
1973  and  operating  level  was  reached  in  August  1974*  Further 
details  of  the  Mica  dam  are  available  in  published  papers  by 
Me i dal  and  Webster  (  1973  )  ,  My Irea  (  1969  ) ,  Webster  (  1  97 0  ), 
and  Webster  and  Lowe  ( 1971),  and  in  unpublished  reports  by 
Caseco  Consultants  Limited* 

At  the  time  of  completion  the  Mica  dam  was  the  highest 
earthfill  dam  in  the  world*  Due  to  its  great  height,  244 
metres,  careful  consideration  was  given  during  design  to 
control  possible  settlements  and  cracking  (Webster,  1970)* 
Extensive  testing  and  analyses  were  carried  out  (Insley  and 
Hillis  1965,  Skermer  1975)  in  order  to  predict  the 
embankment  behaviour  in  advance  of  construction* 

Some  special  features  of  the  design  to  minimize 
cracking  potential  were;  slight  upstream  arching  of  the 
core,  removal  of  river  overburden  material  beneath  the 
central  portion  of  the  embankment,  very  heavy  compaction  of 
the  core  till  and  the  sand  and  gravel  shells*  In  the  event 
of  a  crack  forming  through  the  core  the  excellent 
self-healing  properties  of  the  shell  material  increase  the 
Factor  of  Safety  against  piping  development* 

In  order  that  the  dam  performance  could  be 
monltoredclosely  and  analyzed  an  extensive  instrumentation 
program  was  initiated*  The  instruments  were  Installed  in  4 
main  sections  to  measure  horizontal  and  vertical  movements, 
strains,  earth  pressures  and  pore  water  pressures*  The 
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measuring  devices  included;  82  surface  monuments)  15 
vertical  and  near— vert ica l  movement  gauges,  6  horizontal 
movement  gauges,  207  horizontal  strain  gauges,  53  earth 
pressure  cells,  86  piezometers  and  9  acc e l ogre ph s* 

The  Instrumentation  has  been  carefully  monitored  both 
during  and  subsequent  to  construction*  A  number  of 
back— ana lyse s  have  been  carried  out  since  construction  in  an 
attempt  to  model  the  observed  behaviour  by  finite  element 
methods  (Simmons  1974,  Law  1975,  Skermer  1975)*  In  general 
the  calculated  movements  match  the  observed  movements 
reasonably  well  when  the  correct  boundary  conditions  and 
material  properties  are  entered  in  the  analyses*  No  analysis 
has  been  carried  out  to  predict  the  behaviour  of  the  dam 
during  and  subsequent  to  reservoir  filling* 

2-1,2.  ffiezgBjg.lSj:  Tspeg  ftn<j  L<?Cfttj.Pa6 

A  large  number  of  piezometers  were  installed  in  the 
body  of  the  dam  and  in  the  foundation*  The  foundation 
piezometers  were  placed  both  to  indicate  the  degree  of 
success  of  the  grouting  program  and  to  provide  information 
on  the  seepage  through  the  bedrock  due  to  the  differential 
head  induced  by  the  dam*  The  piezometers  in  the  main  body 
were  installed  to  monitor  the  buildup  of  pore  pressures 
during  construction  and  to  Indicate  the  behaviour  of  the  dam 
during  reservoir  filling  and  during  the  subsequent  years. 

Due  to  the  pervious  nature  of  the  shell  material  only  a 
few  piezometers  (  13)  were  placed  in  strategic  locations  in 
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■these  zones*  to  ensure  that  no  build-up  of  pore  pressure  was 
occurring  in  the  sand  and  gravel  shells*  The  majority  of  the 
piezometers  (32)  were  placed  within  the  till  core*  One  main 
section  and  three  minor  sections  were  instrumented*  The 
sections  are  shown  In  Fig*2*3* 

A  variety  of  piezometer  types  were  installed  in  order 
that  their  reliability  and  behaviour  could  be  compared* 
Having  a  mixture  of  models  could  be  valuable,  in  the  event 
of  one  particular  model  developing  a  problem,  then  at  least 
some  readings  could  be  obtained  from  the  remaining 
piezometers  in  the  section* 

The  types  of  piezometers  installed  were: 

1 •  Standpipe  —  Casagrande  type 
2*  Hydraulic  —  USBR  type 

3*  Pneumatic  —  Hall  type  and  Geosistemas  S« A*  type 
4*  Electric  —  Maihak  type 

These  piezometers  are  described  by  Hanna  (1973),  United 
States  Bureau  of  Reclamation  (1963),  Scott  and  Kilgour 
(1967)  and  their  advantages  and  disadvantages  ere  discussed 
by  Vaughan  ( 1973),  Little  ( 1973)* 

The  choice  of  a  particular  type  of  piezometer  in  each 
location  depended  on  the  following  factors: 

1*  The  dam  zone  or  foundation  material* 

2*  The  distance  from  piezometer  tip  to  readout 
unit  • 

3*  The  relative  elevation  of  piezometer  tip  to 


readout  location* 
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PE  _  Electric  piezometer 

PP  —  Pneumatic  piezometer 

PH  _  Hydraulic  piezometer 
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FIGURE  2  3  Piezometer  Locations 
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Where  long  lead  lengths  were  required  electric  piezometers 
were  installed.  Hydraulic  piezometers  could  only  be  used 
where  the  instrument  house  was  below  the  tip  elevation. 
Standpipe  piezometers  were  located  in  the  cofferdams  and 
were  only  temporary  in  nature.  All  piezomters,  except  the 
standpipes,  were  fitted  with  high  air  entry  porous  stones. 


These  piezometers  were  ordinary  Casagrande  type 
piezometers.  All  these  piezometers  were  abandoned  during 
construction  and  are  not  Important  for  this  thesis.  One  of 
these  standpipes  was  installed  in  the  foundation,  this  was 
destroyed  in  late  1974. 


These  piezometers  are  the  USBR  embankment  type  and  were 
manufactured  by  Advance  Industries  Ltd.,  Vancouver.  They 
were  originally  intended  for  the  Bennett  dam,  B.C. 

13  of  these  were  installed,  fe  in  the  core  and  5  in  the 
downstream  shell.  There  have  been  no  failures  up  to  early 
1979. 

2.  2.3  Engvjtg&lig.  P  iezftmsifcrs 

Two  types  of  pneumatic  piezomters  were  installed.  The 
majority  were  supplied  by  Geo— testing  Inc.,  California  and 
are  otherwise  known  as  Hall  piezometers. 


35  Hall  piezometers  were  installed.  Of  these  6  have 
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failed  and  12  have  been  abandoned  from  the  cofferdams*  The 
remainder  are  still  functioning,  8  in  the  core,  4  in  the 
shells  and  8  In  the  downstream  river  overburden* 

The  second  type  of  pneumatic  piezomter  was  manufactured 
by  Geosistemas  S*A*  of  Mexico*  These  3  piezometers  were 
installed  in  the  upper  core  and  are  still  functioning* 

The  basic  differences  between  the  two  models  are  the 
body  materials  and  the  reading  method*  The  Hall  piezometer 
has  a  stainless  steel  body,  the  Geosistemas  a  PVC  body*  The 
Hall  type  is  read  by  maintaining  a  steady  flow  of  nitrogen 
through  the  system,  whereas  the  Geosistemas  type  is  read  as 
a  pressure  value  upon  closure  of  the  diaphram* 

2 1  2 t  4  Elgclrig  P.i.g£onelen 

The  electric  piezometers  are  the  vibrating  wire  type 
supplied  by  H*Maihak  A*G*  of  Germany*  An  excelllent  review 
of  these  piezometers  is  given  by  Scott  and  Kilgour  (  1967)* 

45  of  these  units  were  installed*  18  have  failed  and  11 
were  abandoned  from  within  the  cofferdams*  Of  the  remaining 
16  one  in  the  core  may  be  suspect*  The  functioning  units  are 
3  in  the  core,  3  in  the  shells,  8  in  the  core  foundations,  1 
in  the  downstream  cofferdam* 

The  location  and  type  of  each  core  piezometer  is  shown 
on  Figure  2*3*  The  vertical  and  near-vertical  movement 
gauges  are  being  utilized  as  open  standpipesas  well  as  their 


normal  function* 
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2.3  Method  of  PI ezomet er 

Two  methods  of  piezometric  installation  are  used  in  the 
Mica  dam,  one  for  the  foundation  piezometers  and  the  other 
for  embankment  piezometers*  Readout  houses  are  located  both 
in  the  upstream  shell  and  the  downstream  shell*  Three  of 
these  houses  located  in  the  upstream  shell  ere  now  flooded 
by  the  reservoir  and  are  abandoned*  One  house  remains  above 
reservoir  level  in  the  upstream  shell*  Three  houses  are 
located  at  various  levels  in  the  downstream  shell  and  it  is 
in  these  that  most  of  the  piezometer  leads  terminate*  Two 
houses  are  located  in  the  abutments*  The  instrument  houses 
are  well  designed  and  are  frost  and  weather-proof t 
containing  ooth  heating  and  lighting* 

The  piezometers  were  read  once  monthly  during  the 
reservoir  impoundment  period.  Subsequent  to  this  they  were 
read  a  minimum  of  3  times  yearly*  The  readings  were 
carefully  tabulated  and  plotted  on  charts  of  level  versus 
time*  The  charts  are  published  In  monthly  or  yearly 
instrumentation  reports  by  Caseco*  Consultants  Limited* 


The  foundation 


are  of  a  limited  interest 


for  this  thesis  and  will  not  be  described  in  detail*  The 


reason  for  this  is  that  too  few  of  the  foundation 
piezometers  are  still  functioning  to  provide  an  accurate  and 
consistent  picture  of  the  exact  conditions  in  the  bedrock* 
These  units  are  simply  installed  in  vertical  drillholes  with 
the  instrument  lines  rising  vertically  out  of  the  drillhole 
and  then  traversing  across  the  embankment*  in  association 
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with  the  other  instrument  leads  in  the  area,  to  the  readout 
house*  The  piezometers  are  sealed  just  above  the  tip  with  a 
bentonite  plug,  and  the  holes  are  carefully  backfilled  to 
ground  level*  The  trench  detail  for  these  instrument  leads 
is  the  same  as  for  the  embankment  type* 

The  embankment  pi ezoae te rs  are  installed  as  the 
embankment  construction  proceeds*  The  tips  are  carefully 
placed  just  below  construction  level*  Then  the  leads  are 
laid  in  a  trench  which  is  offset  from  the  tip*  This  trench 
contains  a  number  of  leads  from  pi ezomete rsa t  different 
locations,  these  leads  are  brought  to  an  instrument  readout 
house*  This  is  normal  procedure  for  embankment  piezometer 
installation  (USER  1963,  Blight  1970). 

The  specifications  for  the  trenches  ore  shown  in  Caseco 
Consultants  Ltd*,  'Contract  documents*  —  vol*2  of  4,  dated 
September  6,  1967*  Details  of  the  layout  are  shown  on  Figure 

2*4*  The  backfill  for  the  core  trenches  was  the  core 
material  with  the  soil  sizes  greater  than  9*5  mm  (3/8  in*  ) 
removed  (Type  A)*  An  alternate  provided  for  was  a  fine 
concrete  aggregate  with  up  to  10%  bentonite  by  volume  added 
(  CSA  A23*  1  )•  The  backfill  was  to  be  hand  compacted  in  75  mm 
(3  in)  layers*  Cut-offs  were  to  be  provided  of  Type  A 
backfill  plus  5%  by  volume  of  bentonite,  not  less  than  0*3  m 
(1  ft)  thick  and  extending  0*3  m  (1  ft)  outside  the  sides  of 
the  trench  at  15  m  ( 50  ft*  )  intervals  or  between  2 
piezometer  tips,  whichever  is  lesser* 

Riser  leads  were  required  for  elevation  changes  between 
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FIGURE  2.4  Instrument-Lead  Trench  and  Riser  Detail 
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lead  trenches.  Details  are  shown  in  Figure  2.7.  Backfill 
should  be  Type  A  material*  with  bentonite  added  every  15  m 
(50  ft.  )  vertically*  to  act  as  a  cutoff. 

The  locations  of  the  piezometer  lead  trenches  and  riser 
pipes  are  shown  in  Figure  2.5.  These  details  were  taken  from 
Caseco  Consultants  Ltd.,  "Dam  Instrumentation  construction 
notes*  1969—1972*  (unpublished).  It  appears  from  the  notes 
that  there  was  some  deviation  from  the  specifications  during 
construction*  which  requires  further  discussion. 

2-sJLtl  Pevjft t  Igo  Enam  &eg.clJE  ica.tigES 

During  1969  the  installation  of  the  leads  in  the 
trenches  closely  followed  the  specifications.  That  is,  5% 
ben toni te— sand  with  a  small  amount  (0.15%)  of  cement  was 
used  150  ram  above  and  below  the  Instrument  leads.  The 
remainder  of  the  trenches  were  filled  with  minus  150mm  till. 
Cut-offs  were  provided  as  specified.  The  risers  were 
backfilled  with  minus  9.5  ram  till. 

During  1970  the  only  reported  change  was  that  the  bulk 
of  the  trench  was  backfilled  with  minus  75  mm  till.  However, 
from  the  limited  data  available  it  appears  that  the  moisture 
content  of  the  backfill  was  very  low,  although  the  densities 
were  reasonably  high. 

Throughout  the  1971  construction  season  the  complete 
trench  was  backfilled  with  sand-ben tonite— cement  mixture 
rather  than  screened  till.  This  was  at  the  contractors* 


choice.  Also  a  statement  was  made  in  the  construction  notes* 
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"it  was  found  that  the  initial  field  moisture 
contents  of  both  Wood  River  gravel  and 
sand-cement-ben  tonite  mix  were  considerably  below 
optimum  and  moisture  was  added  to  improve 
compac  t i on. " 

It  is  apparent  that  some  difficulties  were  encountered  with 

the  risers, 

"Oversize  rocks  and  roots  scalped  from  the  fill  were 
frequently  piled  around  risers  and  the  contractor 
had  to  be  reminded  continuously  to  keep  risers 
clear. " 


The  construction  procedure  for  1972  was  similar  to  that 
of  the  1971  season  discussed  above. 

It  should  be  noted  that  specific  details  of  the  trench 
and  riser  installation  were  difficult  and  in  many  cases 
impossible  to  acquire,  mainly  because  7  years  have  passed 
since  the  construction  was  completed.  However,  the 


"construction  notes"  quoted  from  above  do  provide  some 


insight  into  the  average  method  of  the  installation  and  some 
variations  encountered.  The  significance  of  these  variations 


is  discussed  later  in  Chapter  7 
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CHAPTER  III 


OBSERVED  PIEZOMETRIC  RESPONSE 

Construction  of  the  embankment  commenced  in  April  1969 
and  continued  over  the  next  four  construction  seasons  to 
November  1972  when  the  dam  was  topped  out  at  elevation  763 
m  ,  the  core  depth  being  about  244  m •  Reservoir  filling 
commenced  on  the  29th  March,  1973*  The  "first-filling"  took 
place  over  2  filling  seasons,  this  was  in  order  that  the 
water  level  would  not  be  raised  too  quickly*  A  small 
drawdown  (33  m )  was  observed  between  the  two  filling  stages* 
Elevation  734  m  was  reached  during  this  " f i rs t- f i l 1 1 ng " 
stage  by  mid— August  1974*  Following  this  the  reservoir  has 
fluctuated  seasonally  between  "full  pool"  elevation  (754  to ) 
and  a  winter  level  of  about  730  m*  Maximum  allowable 
drawdown  rates  are  limited  to  about  1*5  m/day  over  a  10  day 
period,  with  much  lower  rates  over  extended  periods  to  a 
minimum  of  0*15  m/day*  Filling  has  averaged  about  0*3  m/day 
over  the  last  5  years* 

The  purpose  of  this  chapter  is  to  discuss  why  the 
piezometers  detailed  In  Chapter  2  have  failed  and  to  briefly 
review  the  response  of  the  piezometers  to  the  reservoir 
filling  and  to  subsequent  reservoir  fluctuations* 
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ILt2  PAg&sg..gtgra  -  EfeULargg  and  Efti lability 

A  considerable  number  of  piezometers  in  the  embankment 


have  either  failed  or  been  abandoned*  The  reason  for  the 
abandonment  is  that  early  in  the  construction  of  the 
cofferdams  some  pizometers  were  installed  to  provide 
additional  design  data*  On  reservoir  filling  many  of  these 
piezometers  have  become  inaccessible,  this  does  not  present 
any  problem  as  the  data  to  be  obtained  was  only  significant 
during  the  early  construction  period* 

The  failures,  however,  were  not  anticipated  and  it  is 
of  interest  to  evaluate  the  reasons  for  these  failures* 

Table  3*1  shows  the  piezometer  category,  number  of  units, 
number  of  failures  and  percentage  success*  It  is  obvious 
that  the  electric  piezometers  have  the  worst  record  by  far* 
From  a  close  inspection  of  the  failure  dates  of  the 
Malhak  units  a  trend  was  immediately  obvious*  Six 
piezometers  failed  in  July  1973*  If  the  problem  was 
straining  of  the  leads,  a  common  cause  suggested  by,  Scott 
and  Kilgour  (1967),  then  the  piezometers  would  not  all  fail 
at  the  same  time*  The  explanation  is  that  two  separate 
lightning  strikes  on  the  dam  caused  voltage  surges  which 
were  too  high  for  the  overvoltage  protection*  This 
phenomenon  has  been  noted  in  the  literature  by  Pinkerton  and 
McDonnell  (1964)  on  the  Tooma  Dam  (Australia)  where  8  out  of 
15  Maihak  piezometers  failed*  Also  Hoskirg  and  Hilton  (1963) 
suggest  that  the  progressive  failure  of  the  electric 
piezometers  in  the  Eucumlene  Dam  were  probably  due  to  a 
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series  of  minor  strikes  because  these  gauges  were  entirely 
unprotected.  It  becomes  obvious  from  a  scan  of  the 
literature  that  the  phenomenon  of  lightning  strikes  is  a 
rare  event t  many  authors  have  had  excellent  success  with 
Maihak  piezometers  and  although  the  literature  does  draw 
attention  to  the  problem  it  does  not  emphasize  lightning  as 
a  potential  danger  (Little  1973,  Scott  and  Kllgour  1967, 
Speedie  1963  ). 

It  is  interesting  to  note  that  after  the  first  failures 
were  diagnosed  in  1970  the  leads  were  altered  in  future 
installation  from  unshielded  ( GG  )  type  to  shielded  G( C  )G 
type.  This  precaution  was  not  sufficient,  however,  because 
of  the  large  number  of  failures  in  1973.  Perhaps  because  of 
the  random  nature  of  lightning  strikes  this  problem  has  not 
been  given  enough  attention,  but  the  consequences  are 
serious  and  more  thought  should  be  paid  to  the  overvoltage 
protect  ion  detail  in  future  inst&llati ons. 

One  piezometer  in  the  core  ( PE35  )  experienced  a  zero 
shift  after  the  lightning  strike  in  1973  and  so  its  readings 
were  re— calibrated  for  this  thesis  after  this  date. 

A  serious  consequence  of  these  failures  is  that  there 
are  virtually  no  reliable  readings  to  be  obtained  from  the 
upstream  side  of  the  core. 

The  pneumatic  piezometers  have  a  moderate  record  and 
the  failures  tended  to  be  progressive  rather  than  sudden. 

The  cause  of  the  failures  has  not  been  established 


definitely  but  are  probably  due  to  excessive  straining  and 
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rupture  of  the  leads  or  leakage  into  the  tubing  over  a  long 
period  (Little  1973,  Vaughan  1973,  Caseco  Consultants  Ltd# 
1970)*  Pneumatic  piezometers  are  a  relatively  new  innovation 
and,  as  such,  their  long  term  reliability  is  not  confirmed 
to  date*  They  may  suffer  from  operater  error,  as  reported  by 
Little  (1973)*  Fortunately,  the  pneumatic  piezometers  are 
not  placed  in  critical  locations  in  the  core* 

The  U*S*B*£2*  type  hydraulic  piezometers  have  an 
excellent  record  in  the  Mica  dam  and  are  considered  to  be 
very  r el i ab l e • Thi s  optimistic  view  of  hydraulic  piezometers 
has  been  shared  by  many  authors  over  the  years* 

The  locations  of  the  functioning  piezometers  are  shown 
on  a  composite  section  in  Fig*3*l* 

All  of  the  piezometers  are  fitted  with  high— ai r-ent ry 
(H*A*E*  )  porous  stones  and  therefore  should  record  the  pore 
water  pressure  in  the  unsaturated  fill*  However,  as  noted  by 
many  authors  (Vaughan  1965,  Scott  and  Kilgour  1967)  in  the 
case  of  an  electric  piezometer  with  no  de-airing  facility 
even  a  H*A*E*  porous  stone  will  not  preclude  air  saturation 
of  the  stone  and  cavity*  In  areas  of  high  pore  pressure  the 
difference  between  the  pore  air  pressure  and  pore  water 
pressure  will  be  small*  However,  in  areas  of  low  pore 
pressures  the  reliability  of  the  electric  piezometers  could 
be  questioned*  Pneumatic  piezometers  also  suffer  from  the 
drawback  of  no  de-airing  facility*  So  the  pneumatic 
piezometers  in  the  upper  core  s  hould  be  treated  with 
caution  as  they  could  very  well  be  reading  pore  air 


Elevation 


25 


780 


E 


740 


Maximum  reservoir 
level  (754m) 


J2. 


/  \ 
PP45  PP46  PP47 


PP50  PP49, PP48 

/•  •/ 

PE44  PP44 


Scale 


0  40  m 

i - 1 - i 


FIGURE  3.1 


Composite 


Section  Showing 


Functioning  Piezometers  (1978) 


26 


pressures*  The  hydraulic  piezometers  have  been  de- aired 
every  few  years  when  necessary  and  so  it  is  felt  that  they 
are  indeed  reflecting  the  pore  water  pressures* 

The  conclusion  is  that  the  hydraulic  piezometers  are 
behaving  excellently  and  the  readings  are  reliable#  The 
electric  units  have  a  very  poor  record  end  the  surviving 
piezometers  may  reflect  pore  air  pressures  in  areas  of  low 
pore  pressure*  One  electric  piezometer  ( PE35  )  is  known  to 
have  suffered  a  zero  shift  in  1973  and  so  is  considered 
unreliable*  The  pneumatic  piezometers  have  a  moderate  record 
and  may  also  be  recording  pore  air  pressures  in  areas  of  low 
pore  pressure*  These  pneumatic  piezomters  are  not  located  in 
critical  areas  of  the  core*  The  nett  result  is  that  the  pore 
pressure  record  for  the  downstream  section  of  the  core  is 
very  good,  whereas  the  upstream  section  is  very  poor* 

InstrMmgflt&tign  Response  In  Mfeigx  ggstlana. 

Dflla.  Saacgg 

The  piezometer  readings  for  the  years  1969  to  1979 
required  for  this  thesis  were  obtained  from  the  Caseco 
Consultants  Ltd*  "Mica  Project:  Reports  on  Dam 
Instrumentation"  1969  —  1979* 

Readings  were  taken  on  an  average  of  once  per  month 
with  the  exception  of  the  winter  when  no  readings  were  made* 
The  monthly  readings  were  taken  from  the  original 
reports  and  combined  in  continuous  plots  over  the  7  year 


period  from  1972  to  1979  (Fig*3*2*  ) 
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□PURE  3  2  General  Piezometer  Levels  versus  Time 
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3*3*2  Ggn&j&A  trend  sl±  rsaefinag 

A  cursory  examination  of  all  the  data  Indicates  a 
similar  pattern  for  all  core  piezometers*  The  upstream  shell 
piezometers  reflect  the  reservoir  head  as  expected  and  the 
downstream  piezometers  show  no  excess  pore  water  pressures* 
The  placement  pore  water  pres  sur  es  we  re  examined*  The 
readings  obtained  within  the  core  were  all  positive  or  very 
slightly  negative*  These  were  all  electric  piezometer 
readings  and  the  reason  could  very  easily  be  that  the 
pizometers  were  reading  the  pore  air  pressure*  At  a 
compaction  water  content  minus  2%  of  optimum  it  is  expected 
that  fairly  large  negative  pore  presures  would  develop*  It 
has  been  noted  previously  by  Little  (1973)  that  electric 
piezometers  may  not  reflect  this  negative  pore  water 
pressure  successfully*  The  reason  that  no  hydraulic 
piezometric  measurements  were  available  immediately  after 
placement  was  that  the  instrument  house  installations  were 
not  prepared  at  that  time* 

Next  the  end  of  construction  pore  pressures  were 
inspected*  As  expected  with  the  increased  total  stresses  due 
to  the  application  of  the  remaining  fill  the  pore  pressures 
had  become  positive*  Contours  of  the  end  of  construction 
pore  pressures  are  shown  in  Fig*3*7* 

Some  consolidation  is  apparent  between  the 
end— of-cons truct i on  and  first  filling  in  all  the  lower 
piezometers*  The  picture  is  confused  in  the  upper  elevations 
because  the  applied  overburden  stresses  were  not  large 
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enough  to  indicate  specific  pore  pressure  gradients,  coupled 
with  probable  ambiguity  over  pore  air  and  pore  water 
pressures#  (Figure  3*3)  One  unusual  feature  noted  during 
this  phase  is  that  the  central  core  zones  appear  to 
consolidate  more  rapidly  than  the  outer  core  zones  which 
might  be  unexpected* 

"F irs t— f i ll i ng"  of  the  reservoir  was  reflected  in  the 
lower  core  piezometers  by  a  small  increase  in  pore 
pressures*  However,  during  the  second  stage  of  the 
"first— filling"  there  is  a  very  noticeable  increase  in  the 
pore  pressure  response*  This  change  in  response  occurred 
approximately  1  1/4  years  after  filling  commenced*  Following 

the  next  filling  peak  some  of  the  piezometers  show  a  d 
pressure* 

From  September  1974  onwards  is  termed  the  'steady 
reservoir  level'  even  though  there  are  seasonal  fluctuations 
of  approximately  24  m*  The  lower  core  piezometers  tended  to 
model  these  fluctuations  with  a  time  lag,  the  upper 
piezometers  also  show  some  response  to  the  reservoir 
f luctuat i ons • 

The  general  conclusion  was  that  perhaps  the  piezometers 
were  reacting  to  the  two  processes  involved  during  reservoir 
filling  and  transient  flow  leading  to  "steady  state" 
conditions*  These  two  processes  are: 

1*  Total  stress  change* 

2*  Transient  seepage* 

The  term  "steady  state"  refers  to  a  condition  which  is  close 


, 
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FIGURE  3.3  Contours  of  End -of-Construction  Pore  Pressures 


31 


to  the  long  term  steady  state*  This  state  may  never  be 
completely  reached  while  the  resevoir  fluctuates* 


.3 1.4  Ingtrygeal^tlop  Be.talla  and  Behaviour  Ellhlp  Ma.iar 
Sections 


3t  Sepllon  22+50 

This  is  the  major  instrumentation  section*  The 
of  the  piezometers  is  shown  in  Fig*2*3*  There  are  4 
piezometers  at  elevations  590*7  m»  650*7  m,  706  in,  a 
m.  The  upper  two  rows  contain  2  piezometers  each,  an 
lower  two  rows  contain  4  piezometers  each*  Theoretic 
this  would  have  supplied  an  excellent  indication  of 
pressure  distribution  within  this  section*  Unfortuna 
number  of  piezometers  either  failed  or  their  reading 
invalidated  for  other  reasons  discussed  previously* 

Piezometer  PE23  failed  in  August  1970,  PE25  and 
failed  in  July  1973,  the  remainder  performed  satisfa 
However,  the  readings  from  PE38,  PP36  and  PE44  are 
invalidated  because  the  leads  from  these  piezometers 
Installed  in  trenches  which  ran  north  into  the  upstr 
shell  and  these  piezometers  are  essentially  reflectl 
upstream  water  levels*  From  a  total  number  of  12  pie 
only  5  are  considered  to  be  functioning  and  reliable 
these  piezometers  are  located  in  the  downstream  sect 
the  dam  core  with  the  exception  of  PE24*  Water  level 
beenread  continuously  in  the  vertical  movement  gauge 
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MV8,  MV9*  These  movement  gauges  have  been  founded  in  the 
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bedrock,  and  they  permit  the  entrance  of  seepage  along  their 
length*  For  this  reason  they  should  reflect  the  phreatic 
surface  location  with  time*  It  is  unlikely  that  high  bedrock 
seepage  pressures  are  distorting  the  readings  because  the 
gauges  are  grouted  into  the  bedrock* 

E,\eYatl<2S  m 

Piezometer  tips  PE24,  PHI  and  PH2  are  located  at  this 
elevation*  These  three  piezometers  show  a  similar  trend, 
each  reflecting  the  consolidation  in  the  core  until  the  pool 
elevation  raises  above  590*7  m  when  a  continuous  but  slight 
increase  in  pore  pressure  is  shown*  During  the  first 
drawdown  in  the  winter  of  1973,  PE  24  and  PHI  continue  to 
rise  whereas  PH2  decreases  slightly*  During  the  second 
filling  the  3  piezometers  respond  more  strongly  and  from 
this  time  on  they  all  respond  in  a  similar  but  damped  manner 
to  the  reservoir*  At  all  times  the  pore  pressure  at  PE24  is 
greater  than  PHl  which  is  in  turn  greater  than  PH2*  In 
principal  the  trend  is  correct  and  the  magnitudes  will  be 
investigated  in  following  sections* 

Elevation  a 

Piezometer  tips  PH6,  PH7  are  installed  at  this 
elevation*  The  first  notable  characteristic  of  these  two 
piezometers  is  that  their  respective  pore  pressures  are  very 
close  in  magnitude  and  they  fluctuate  in  a  very  similar 


manner* 


, 
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They  both  continue  to  decrease  at  an  average  rate  of 
0.5  m/mo.  until  the  second  filling  stage  when  they  begin  to 
react  to  the  reservoir  level.  They  continue  to  rise  until 
the  'steady*  reservoir  condition  is  reached  in  1976  when 
they  appear  to  stabilize  and  follow  the  reservoir  level.  The 
difference  in  readings  between  PH6  and  PH7  appears  to  be  too 
small  to  reflect  normal  behaviour. 


3  •  .4  .t  2 


24+5Q 


This  is  the  only  other  Instrumented  deep  section.  It 
contains  3  piezometers  and  these  were  all  installed  close  to 
the  bedrock/fill  interface  (see  Fig. 2. 3).  One  of  these, 

PE22,  failed  in  August  1970.  The  remaining  two  PP20  and  PP21 
functioned  until  July  and  August  1974,  respectively.  It 
would  appear  that  these  piezometers  were  reflecting  rather 
high  pore  pressures  and  this  may  have  been  due  to  the 
proximity  of  the  foundation  boundary.  Their  more  rapid  and 
more  extreme  response  at  the  early  stage  of  filling  Is 
probably  due  to  the  higher  degree  of  saturation  of  the  core 
material  at  this  level,  which  is  a  function  of  the  higher 
stresses,  seepage  from  bedrock,  and  seepage  from  riverbed. 

The  most  likely  cause  of  failure  was  the  rupture  of  the 
leads  from  the  pneumatic  piezometers,  caused  by  excessive 
straining  of  the  tubes. 


3*4.3  Sec  tlQii  lfr+5Q 


This  section  contains  2  rows  of  piezometers, 
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PE35***PH10  at  elevation  677  m  and  PE36***PH11  at  elevation 
716*6  in*  PE40  (which  is  located  below  PH10)  and  PE36  both 
failed  in  July  1973* 

From  a  total  of  5  piezometers,  2  have  failed*  This 
section  is  important  though  because  PE35  is  still 
functioning  and  it  is  located  close  to  the  upstream  edge  of 
the  core* 


Piezometers  PE35  and  PH10  are  located  at  this 
elevation*  During  July /August  1973  PE35  suddenly  dropped 
about  6*5  m*  This  drop  was  not  reflected  by  any  other 
piezometer  and  the  reservoir  was  still  rising*  It  is 
apparent  that  this  drop  of  6*5  m  is  a  direct  consequence  of 
the  electric  storm  which  caused  the  other  electric 
piezomters  to  malfunction*  Another  observation  to  indicate 
the  malfunction  is  that  piezometer  PH10  which  is  downstream 
of  PE35  was  actually  reading  higher  than  PE35  during  1975* 
However,  the  general  trend  of  PE35  seems  correct  during  the 
following  years,  and  it  is  reasonable  to  assume  that  the 
lightning  strike  may  have  only  caused  a  zero  shift  in  the 
readings*  So  the  readings  have  been  corrected  upward  by  6.5 
m*  Any  conclusions  drawn  solely  from  PE35  though,  should  be 
treated  with  due  caution* 

PH10  performs  normally*  The  only  unusual  occurrence  is 
a  sudden  drop  of  approximately  0*75  m  on  Dec* 1  1976*  This 


was  caused  by  a  gauge  recalibration  at  that  time*  Otherwise 
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PH10  is  quite  similar  in  form  to  the  other  piezometers# 

Elgyallgn  ?16jfr  m 

PHll  is  similar  in  most  respects  to  PH10#  Including  a 
re cal i bra ti on  on  December  1,  1976# 

3. 4#  4  Sec tian  30+50 

This  section  contains  2  rows  of  piezometers#  PE34t  PH9 
at  elevation  677  m  and  PE12  at  elevation  718  m#  Included  in 
this  section  is  PE39  which  is  projected  from  Section  31+15# 
PE34»  PE37  and  PE39  all  failed  in  August  1973  and  once  again 
all  these  piezometers  were  located  on  the  upstream  side  of 
the  core# 

Elevation  677  m 

PH9  is  functioning  normally#  It  does  not  begin  to  react 
until  the  second  filling  phase  in  1974  and  then  it  slowly 
begins  to  model  the  reservoir  behaviour# 

Eleva  tion  718  El 

PH 12  also  appears  to  be  functioning  normally#  It  reacts 
only  very  slightly  until  mid— 1976  when  it  too  models  the 
reservoir  behaviour#  This  is  expected  because  of  the 
relatively  high  tip  elevation# 

Some  individual  piezometers  are  located  between 
sections  14+00  to  15+30  and  34+35  to  36+05#  These  are 
piezometers  PP48fPP49f  PP50  and  PP45,  PP46,  PE45,  PP47, 
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respectively*  The  tip  elevations  of  these  groups  are  at 
approximately  747  m  and  750*7  m* 

The  reservoir  level  does  not  reach  elevation  747  m 
until  mid— 1976;  it  remains  above  this  level  for  6  months*  it 
does  not  cross  this  elevation  again  until  mid-1978,  except 
for  a  month  in  1977  when  the  reservoir  peaks  at  748  m* 

As  expected*  these  piezometers  show  very  limited 
response,  the  upstream  piezometer  PP50  showing  more 
fluctuations  than  the  others* 

These  piezometers  have  not  been  subject  to  detailed 
analysi s* 


Co, 


The  installation  of  32  piezometers  in  the  embankment 
core  has  resulted  in  12  failures,  4  completely  unreliable 
( PE38 * PE36 * PE44, PP44  )  *  1  may  be  unreliable,  leaving  a 

remainder  of  15  piezometers  functioning  normally*  The 
distribution  of  the  reliable  piezometers  1e  unfortunate  in 
that  only  2  of  these  reflect  the  pore  pressure  behaviour  in 
the  upstream  side  of  the  core* 

The  general  behaviour  of  the  piezometers  appears 
normal;  however,  this  chapter  did  not  investigate  the 
magnitudes  of  the  pore  pressures  which  are  dealt  with  in 
later  chapters  in  more  detail* 

The  very  poor  performance  with  electric  piezometers 
could  probably  have  been  avoided  by  the  installation  of  more 
reliable  overvoltage  protection  devices*  The  installation 


, 
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method  Itself  apppears  to  have  been  goo d  because  there  are 
no  failures  due  to  excessive  straining  or  rupture  of  the 
leadst  a  common  fault  noted  by  Scott  and  Kilgour  (1967). 

The  performance  of  the  pneumatic  piezometers  has  not 
really  been  conclusive  in  the  embankment*  becuase  the 
majority  of  these  piezometers  are  at  very  high  elevations 
where  the  responses  are  not  easy  to  analyze. 

The  hydraulic  piezometers  have  behaved  excellently.  The 
reliability  of  hydraulic  piezometers  has  been  noted  by  many 
authors  (Little  1973*  Vaughan  1973,  Mackellar  et  al  1973* 
Hosking  and  Hilton  1963)  provided  the  installation  is 
carefully  handled  by  experienced  personnel. 

The  reccommendat i on  based  on  the  Mica  Dam  performance 
is  that  in  critical  locations*  where  possible,  hydraulic 
piezometers  should  be  used. 
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CHAPTER  IV 


STEADY  STATE  SEEPAGE  ANALYSIS 
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One  of  the  major  concerns  of  this  thesis  is  to  study 
the  process  leading  to  steady  state  flow  conditions  in  the 
Mica  Dam*  In  order  to  ascertain  this  it  is  necessary  to  know 
what  the  theoretical  steady  state  conditions  are  for  this 
dam*  A  number  of  methods  are  available  to  analyse  the  dam 
for  steady  state  conditions* 

These  are  as  follows? 

1*  Closed  form  solutions 

2*  Electric  analogue  or  viscous  fluid  models 
3*  Graphical  solutions  (  flow  nets) 

4*  Approximate  solutions  by  numerical  methods 
( finite  difference  or  finite  element) 

Graphical  solutions  are  probably  the  simplest,  but  because 
of  the  large  number  of  flow  nets  which  would  be  required  and 
the  difficulties  encountered  with  the  seepage  face,  this 
method  was  not  adopted*  The  most  suitable  method  appeared  to 
be  by  finite  element  techniques*  A  computer  program  already 
existed  within  the  department  which  could  handle  this 
problem*  This  program  (FPM5)  was  developed  by  Taylor  and 
Brown  at  the  University  of  California  in  1967* 
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4.J.2.  El.nl.ta  Element  £aa jgular  Program  r  EEM5 

The  finite  element  seepage  program  used  during  this 
thesis  was  originally  proposed  by  Taylor  and  Brown  (  1967  )# 
The  program  has  been  used  previously  at  the  University  of 
Alberta  by  Guther  (1972b)  during  his  Master  of  Science 
thesis  work#  Guther  (1972a)  prepared  the  Users  Manual,  Soil 
Mechanics  No#  16  entitled  " Axi syrame tr ic  and  Plane  Flow  in 
Porous  Media",  which  describes  the  preparation  of  input  and 
use  of  the  program#  Further  details  on  the  use  of  this 
program  are  given  by  Kealy  and  Busch  (  1971  )• 

The  finite  element  theory  for  use  in  seepage  problems 
through  porous  media  has  been  adequately  described  elsewhere 
(Zienkiewicz  1971,  Taylor  and  Brown  1967,  Guther  1972),  and 
so  will  not  be  further  discussed  here# 

The  program  FPM5  deals  with  the  flow  of  liquids  through 
a  saturated  porous  medium  where  motion  is  governed  by  a 
generalized  Darcy  equation#  It  has  the  capability  of  dealing 
with  a  free  surfaceand  ax isymmet ric  or  plane  flow  problems# 
The  finite  element  method  is  very  versatile  and  can 
easily  handle  effects  such  as  anisotropy,  non-homogeneity 
and  difficult  boundary  geometries#  The  accuracy  of  the 
program  is  indicated  by  Guther  (1972)  in  a  comparison  with  a 
closed  form  solution# 

An  ambiguity  does  exist  at  the  point  of  intersection  of 
the  free  surface  with  the  boundary,  because  the  nodal  flow 
at  this  point  cannot  be  zero#  It  is  important  to  reduce  the 
mesh  size  in  this  area  to  minimize  this  ambiguity  effect. 
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and  so  reduce  the  error*  This  problem  is  discussed  by  Neuman 
and  Witherspoon  (  1970)  and  they  show  a  case  where  the 
solution  diverges  rather  than  converges;  however*  it  is  felt 
that  with  careful  attention  to  the  mesh  size  in  this  area 
the  program  would  then  tend  to  converge  on  the  correct 
answer* 

The  core  in  the  Mica  dam  overhangs  the  downstream 
shell*  This  results  in  an  unsaturated  zone  below  the  seepage 
face  in  the  downstream  shell  (Casagrande,  1937)*  The  FPM5 
program  cannot  deal  with  this  situation  so  the  core  has  been 
treated  as  a  single  unit  with  appropriate  boundary 
condl t ions • 


±a.. 2  Ennndgca:  £gniHtlanj5  and  Ma.tg.rlal  Exmngrllg.fi  a  £  i.hs.  Mlaa 
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A  general  section  of  the  Mica  Dam  is  shown  in  Fig*2*2* 
The  material  properties  of  the  different  zones  are  discussed 
in  detail  by  Law  (1975)*  Simmons  (1974),  Meldal  and  Webster 
(1973),  and  Webster  (1970)*  Some  till  properties  are  listed 


below : 

Liquid  limit  ••••••  18% 

Plastic  limit  •••••  13% 


Plasticity  index  •••  5% 

The  most  important  parameters  required  for  the  seepage 
analysis  were  the  relative  permeabilities  of  the  different 
zones,  the  effectiveness  of  the  foundation  grouting,  and  the 
degree  of  anistropy  within  the  compacted  fill*  These 
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properties  have  been  reported  in  Caseco  Reports  (various), 
Simmons  (  1974), Law  (  1975  ),  with  the  exception  of  the  degree 
of  anistropy  o-f  the  fill. 

The  shell  material  has  a  permeability  of  about 
10— 3cm/sec.  which  is  approximately  3  to  4  orders  of 
magnitude  larger  than  the  core  permeability.  It  is 
reasonable  then  to  treat  the  core  as  a  separate  unit  because 
the  shells  are  essentially  ‘free-draining*  and  will  not  have 
any  influence  on  the  pore  pressure  distribution  within  the 
core  . 

The  permeability  of  the  core  fill  is  a  function  of  many 
variables  such  as  particle  size*  void  ratio,  composition, 
fabric  and  degree  of  saturation  (  Lambe  and  Whitman,  1969) 
or,  in  other  terms,  the  permeability  is  a  function  of  the 
soil  structure  which  in  turn  is  influenced  by  the  compaction 
method,  compaction  water  content  and  effective  confining 
pressure.  These  values  were  obtained  from  direct 
permeability  tests  during  extensive  triaxlal  testing  on 
representative  samples  reported  by  Insley  and  Hillis  (1965). 
A  figure  (Law  1975)  has  been  extrapolated  to  provide  values 
of  permeability  At  the  low  compaction  water  contents 
observed  in  the  upper  core.  These  results  may  be  slightly  in 
error  but  they  do  indicate  a  lower  bound,  see  Fig.  4.1. 

The  minor  principal  effective  stress  is  estimated  to 
range  from  500  kPa  in  the  lower  core  to  350  kPa  in  the  upper 
core  (Simmons  1974).  Using  these  values  it  Is  assumed  that 
the  lower  core  has  an  average  permeability  of  10— 7cm/sec. 


- 


COEFFICIENT  OF  PER MEA  Bl  LI  T  Y  —  cm/sec  x  10~ 


FIGURE  4.1  Core  Permeability  vs.  Compaction  Moisture  Content 

(after  Law  1975) 
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and  th«  upper  till  at  the  lower  water  content  has  a 
permeability  of  10-6cm/sec*  or  greater*  The  reason  for  this 
extrapolation  is  justified  by  examination  of  the  figures  in 
the  paper  by  Mitchell  et  al  ( 1965)  which  clearly  shows  the 
influence  on  permeability  of  a  wider  range  of  moisture 
contents*  Webster  ( 1970)  indicates  the  till  permeability  to 
be  about  10— 7cm/sec*  and  the  sand  and  gravel  shells  to  be 
between  2  x  10— 3cm/sec*  and  9  x  10— 3cm/sec*  depending  on  the 
dry  densities* 

The  degree  of  anistropy  of  the  till  core  has  never  been 
evaluated  but  It  is  unlikely  to  be  greater  than  4  (Sherard 
et  al*  1963)*  Values  as  high  as  10  were  used  in  this 
analysis,  which  is  certainly  an  upper  bound* 

The  seepage  through  the  bedrock  foundation  was  more 
difficult  to  evaluate*  Three  piezometer  sections  were 
installed  to  obtain  information  on  the  effectiveness  of  the 
grouting  program*  Unfortunately,  a  number  of  the  piezometers 
have  failed  and  the  information  from  the  remaining 
piezometers  is  rather  confusing*  It  would  appear  that  the 
grouting  program  was  not  wholly  successful  in  reducing  the 
underseepage*  Sensitivity  of  the  pore  pressure  distribution 
to  this  boundary  condition  was  investigated* 

The  assumption  with  respect  to  the  shells  was  much  more 
s  t  raj.  ght  forward*  The  piezometers  located  in  the  upstream 
shell  indicate  that  the  head  loss  in  this  zone  is  indeed 
negligible*  Therefore,  the  reservoir  level  is  assumed  to  act 


directly  on  the  upstream  face  of  the  core*  The  piezometers 
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in  the  dowstream  shell  indicate  that  this  material  is 
1  free-draining'  »  seme  small  positive  pore  pressures  are 
probably  due  to  surface  run-off  and  infiltration*  The  lower 
piezometers  give  a  reasonable  Indication  of  the  •tailwoter* 
elevation  in  this  shell* 

Figure  4*2  shows  the  model  assumed  for  the  steady  state 
seepage  analysis* 

kgc.ft t i<?n  o_f  the  downstream  water  leye  1 
The  flow  in  the  downstream  shell  was  assumed  to  be 
horizontal  or  sub— ho r i zo nta l »  so  the  pore  pressures  in  the 
piezometers  should  indicate  the  phreatic  surface  in  this 
zone*  For  the  main  Section  22+50  the  readings  inaicate  e. 
phreatic  surface  which  is  more  or  less  constant 

(  fluctuations  less  than  6  m)  the  fluctuations  probably  being 
due  to  rainfall  and  run-off*  There  is  a  consistent  drop  in 
phreatic  surface  towards  the  toe  and  the  lower  piezometer 
PP23  has  a  greater  pore  pressure  than  PP26  which  may 
indicate  a  slight  upward  flow  from  the  bedrock • Thes e 
variations  are  not  significant* 

The  phreatic  surface  in  the  downstrean  shell  is  taken 
at  573  m  (+  or  3  m)  with  some  confidence*  It  may  be  seen 
later  from  the  analysis  that  the  pore  pressure  distribution 
is  not  very  sensitive  to  this  exact  elevation* 

Sections  16+50  and  30+50  are  less  precise*  There  are  no 
piezometers  located  in  the  downstream  shell  below  the 
phreatic  surface*  PH14  and  PH15  are  both  above  this  level* 
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CORE  MODEL 


FIGURE  4.2  Model  for  Steady  State  Seepage  Analysis 
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The  phreatic  surface 
which  results  from  a 
nature  of  the  valley 


assumed  (65*6  m)  is  relatively 
consideration  of  the  3  dimensie 
and  the  high  shell  permeability 


oJ.  Ihs. 


W&l££  l£Y£.l 


As  discussed  previously  this  level  is  assumed 
equivalent  to  the  reservoir  level*  The  piezometers  i 
zone  substantiate  this  assumption*  as  do  the  vertica 
movement  gauges*  So  there  is  negligible  head  loss  ac 


upstream  shell  and  there  is  no  time  lag  involved 


dLt.3-«.3  g.gg.pQfi.e  Iknoygh  lh£.  t>£d£ack  £gppdfrti gp 

The  bedrock  beneath  the  core  is  described  as  "a 
metamorphic  series  of  late  Precambrian  age  with  main 
types  of  mica  schists  and  granite  gneisses"  (Meidal 
Webster  1973)*  Joint  spacing  is  often  less  than  0*6 
the  surface  increasing  with  depth  to  more  than  15  m* 
testing  in  the  exploration  drill  holes  indicated  the 
be  "relatively  tight"  except  near  the  surface* 

A  grouting  program  was  considered  necessary*  Bl 
grouting  was  carried  out  to  a  depth  of  9  m  on  a  6  m 
grid*  In  addition*  a  line  of  deep  (approx*  30  m)  hoi 
apart  were  grouted*  The  effectiveness  of  these  roeasu 
monitored  by  piezometers  (Fig*4*3)* 

Section  22+50  indicates  grouting  to  be  about  30 
effective*  Section  31+15  indicates  about  85%  effect! 
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Section  16+50  showed  the  grouting  to  be  ineffective  up  to 
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Scale  l. 


40  m 
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FIGURE  4.3  Bedrock  Model  and  Piezometer  Locations 
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the  piezometer  failures  in  July  1973* 

The  Instrumented  sections  indicate  the  possibility  of 
significant  underseepage*  So  sensitivity  to  this  boundary 
condition  was  investigated* 

A  permeability  value  for  the  bedrock  was  assumed  of  the 
order  of  10— 3cm/sec*,  there  is  no  accurate  basis  for  this 
assumption  but  it  is  considered  to  be  an  upper  bound  value* 


±iA  Eulalia  and 


Ql  Analysis 


Finite  element  meshes  were  drawn  for  the  four  main 
instrument  sections*  The  meshes  were  varied  for  different 
reservoir  levels*  Initially,  the  piezometer  locations  were 
not  modelled  by  fixed  nodes*  However,  to  increase  the 
accuracy  it  was  later  decided  to  provide  fixed  nodes  at  each 
piezometer  location  where  practical*  It  should  be  noted  that 
in  the  early  stages  of  analysis  PE24  was  located 
incorrectly,  too  close  to  PHI*  Because  of  an  erroneous 
statement  in  a  Caseco  Instrumentation  report  dated  1S70, 
"•••pore  pressure  values  in  PHl ,  which  is  installed 
very  close  to  PE24,***" 

This  report  is  confusing  PE23  with  PE24*  The  location  error 
was  observed  during  later  inspection  of  construction  notes, 
and  the  location  of  PE24  was  corrected* 

The  value  of  permeability  chosen  initially  for  the 
analyses  was  10— 6cm/sec*  The  absolute  magnitude  of 


permeability  does  not  have  to  be  exact  for  the  analysis. 


unless  the  flow  values  are  of  interest*  It  is  the  various 
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permeability  ratios  which  control 
within  the  core*  The  permeability 
analysis  are  interpreted  from  the 
(  Law,  1  976  )• 


the  pore  pressure  state 
ratios  used  in  the 
published  labatory  testing 


4j.4_?.l  Ssclign  22±fL£ 

Firstly,  the  sensitivity  of  the  analysis  to  the 
downstream  water  level  was  checked*  A  range  of  levels  from 
527  m  to  597  m  was  chosen  for  a  constant  reservoir  level  at 
734  m*  In  the  range  of  interest,  573  +  or  -3m,  the 
variation  in  piezometric  levels  at  elevation  590*7  ra  (PE24, 
PHI,  PH2 )  was  less  than  0*5  m*  At  the  higher  elevations 
there  was  no  significant  variation*  The  conclusion  was  that 
moderate  variations  in  downstream  water  elevation  have  only 
a  local  effect  on  the  pore  pressure  distribution  in  the 
vicinity  of  the  water  elevation* 

Next  the  accuracy  of  the  solution  was  checked  with 
respect  to  mesh  size  and  tolerance  value*  A  fine  mesh  (184 
elements)  and  a  coarse  mesh  (120  elements)  were  compared* 
The  coarse  mesh  provided  results  within  2%  of  the  finer 
mesh*  This  was  considered  accurate  enough  because  a  slight 
change  in  tip  location  would  cause  an  error  of  this 
magnitude*  Of  course  using  the  coarse  mesh  caused  the 
boundary  ambiguity  effect  to  be  more  pronounced* 

The  tolerance  level  was  checked  on  the  fine  mesh  by 
altering  the  tolerance  from  0*2  m  to  0*02  m*  This  increased 
the  number  of  iterations  required  from  6  to  15  with  no 
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significant  change  in  phreatic  surface  level#  The  results  of 
these  analyses  are  shown  on  Fig*4*4*  The  conclusion  was  that 
the  coarse  mesh  is  satisfactory  and  a  tolerance  level  of 
0*1%  of  the  available  head  is  adequate  to  give  results  of 
acceptable  accuracy* 

A  number  of  computer  runs  were  made*  with  the 
assumptions  discussed  above f  for  varying  reservoir  levels. 
The  core  was  considered  homogeneous  and  isotropic*  A  sample 
of  the  resulting  grids  and  equipotent ial  distributions  are 
shown  in  Figs*4*5  to  4*7*  The  total  head  was  plotted  versus 
reservoir  level  for  each  piezometer*  The  results  are  shown 
on  Fig*4*8* 

Next  the  effect  of  anistropy  was  investigated.  A 
horizontal  to  vertical  permeability  ratio  of  10  was  chosen 
as  an  upper  bound*  The  effect  of  this  anistropy  was  very 
slight  indeed*  The  phreatic  surface  was  close  to  horizontal 
but  the  piezometric  readings  were  not  effected  in  the  lower 
areas*  The  area  of  major  difference  is  the  area  close  to  the 
phreatic  surface  in  the  downstream  section*  In  other  words » 
the  anisotropy  causes  only  a  local  disturbance  of  pore 
pressure  distribution*  Sample  results  are  shown  in  Fig*4*9* 

The  effect  of  non-homogeneity  was  investigated  for  the 
differing  core  properties  at  elevation  610  m*  Above 
elevation  610  m  the  permeability  chosen  was  10-6cm/sec* 
below  this  elevation  10— 7cm/sec*  This  effect  was  also 
studied  for  isotropy  and  anisotropy*  When  the  reservoir  is 
at  higher  elevations  (  734  m  )  the  disturbance  is  not 
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FIGURE  4.4  Sensitivity  of  Analysis  to  Mesh  Size  and  Tolerence 
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FIGURE  4.5 


Steady  State  Analysis.  Section  22»50.  Res.  Level  754m. 
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FIGURE  4.6 


Steady  State  Analysis.  Section  22 ♦  50.  Res.  Level  734  m. 
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EIGURE  4.7  Steady  State  Analysis.  Section  22*50 .  Res.  Level  660m. 
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FIGURE  4.8  Calculated 


Steady  State  Piezometer  Readin 
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FIGURE  4.9 


Effect  of  Material  Anisotropy.  Section  22*50 


57 


pronounced  but  remains  quite  local  in  the  area  of  the 
material  change*  However,  for  lower  reservoir  levels  (  660  ®  ) 
the  alteration  in  pore  pressure  distribution  is  more  severe 
but  still  only  locally*  The  net  result  is  that  the 
plezometric  readings  above  elevation  590  m  do  not  exhibit 
any  noticable  change  in  readings  due  to  this 
non-homogeneity*  The  results  are  shown  on  Fig*4*10* 

The  final  boundary  condition  investigated  was  concerned 
with  the  foundation  bedrock  seepage*  A  simple  model  was 
chosen  with  reservoir  level  at  660  m*  Isotropy  was  assumed* 
Three  material  types  were  chosen: 

1*  core  material:  k  =  10-6cra/sec* 

2*  Bedrock  :  k  =  10— 3cm/sec* 

3*  Grouted  section  :  k  =  10— 4  cm/sec* 

This  model  assumes  the  grout  curtain  to  be  efficient  but  the 
bedrock  below  to  be  of  high  permeability* 

The  resultSf  Flg*4*ll  compared  to  Fig*4*7»  indicate  a 
local  change  in  potential  distribution  around  the  base  of 
the  core,  close  to  the  grouted  section*  No  effect  is  felt  at 
higher  elevations  in  the  core*  Although  this  model  is  very 
crude  and  the  behaviour  is  not  closely  reflected  by  the 
foundation  piezometers  it  was  felt  that  since  this 
underseepage  has  only  a  localized  effect  on  the  core  pore 
pressure  distribution,  it  was  not  worthwhile  refining  the 
model*  The  results  are  compared  with  the  piezometric 
readings  in  Flgs*4*12* 


The  conclusion  from  this  parametric  study  is  that  none 
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FIGURE  4,10 


Effect  of  Material  Non -Homogeneity.  Section  22*50 
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FIGURE  4.11 
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Notes 


Scale : 


0 

I_ L 


1.  Reservoir  elevation  660m  for  model. 

2.  Downstream  water  level  573m  for  model. 

3.  All  piezometers  referred  to  grout  curtain  as  origin. 

4.  Figures  in  brackets  (  )  are  measured  potentials. 

5.  Measured  potentials  are  for  a  reservoir  level  at  732m. 


40  m 


FIGURE  4.12  Comparison  of  Bedrock  Model  and  Piezometric  Readings 
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of  the  conditions  investigated  change  the  pore  pressure 
distribution  within  the  core  in  a  significant  manner,  with 
the  exception  of  very  local  effects*  The  use  of  hoegeneous 
isotropic  material  properties  will  be  sufficiently  accurate* 

4*4x3  S&£llon  24±5£ 

Very  little  analysis  was  done  on  this  section  because 
of  the  proximity  of  the  piezometers  to  the  bedrock*  hence* 
the  localized  effects  of  underseepage  became  important* 
Results  of  an  analysis  are  shown  on  Fig«4«13«  But  they  are 
probably  incorrect  near  the  base  of  the  core  due  to  the 
uncertain  boundary  condition* 

4«  4f  3  tlgfl  16±5£ 

This  section  was  analyzed  as  a  homogeneous  isotropic 
section*  for  various  reservoir  levels*  The  section  is 
entirely  above  elevation  6T 0  m*  So  there  is  ro  abrupt  change 
in  core  properties*  The  results  are  shown  on  Fig*4*l4*  And  a 
plot  of  piezometric  level  versus  elevation  was  made 
(Fig*4*8)*  An  anisotropic  core  was  analyzed  to  investigate 
the  local  effect  on  piezometer  PHll*  The  result  is  shown  in 
Fig*4*15*  There  is  a  moderate  local  disturbance,  the 
piezometric  level  is  raised  by  1*8  m* 

4*4*4  SeQr.tlgn  2L£±5i2 

A  minimum  number  of  cases  were  studied,  due  to  the 
similarity  between  this  section  and  Section  16+50*  The 
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FIGURE  413 


Steady  State  Analysis.  Section  24*50 
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FIOURE  4.14  Steady  State  Analysis.  Section  16*50 
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Effect  of  Material  Anisotropy.  Section  16*50 
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results  and  conclusions  were  very  similar*  Results  are  shown 
on  Figs*  4*16* 


This  analysis  has  enabled  plots  of  reservoir  level 
versus  piezometer  level  to  be  drawn  for  the  major 
instrumentation  sections  for  steady  state  seepage 
conditions*  Various  effects  such  as  enisotropy, 
non-homogeneity,  underseepage  were  investigated  and  found 
only  to  have  localized  effects  which  although  generally 
unimportant  should  not  be  ignored  in  certain  cases,  for 
example,  when  piezometers  are  located  close  to  geometric 
boundaries  with  ill— defined  boundary  conditions*  In  order  to 
increase  the  accuracy  of  the  analysis  it  would  be  necessary 
to  have  many  more  specific  details  on  such  parameters  as 
construction  irregularities  and  anisotropy  which  are  not 
readily  available*  However,  it  is  considered  that  the 
results  are  well  within  an  acceptable  degree  of  accuracy  and 
should,  indeed,  compare  with  actual  piezometer  readings  when 
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FIGURE  4.16  Steady  State  Analysis.  Section  30+50 
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CHAPTER  V 


COMPARISON  OF  OBSERVED  AND  PREDICTED  PORE  WATER  PRESSURES 

5jl1  J  Plug  due  tlQH 

The  purpose  of  this  chapter  is  to  compare  the  actual 
pore  pressure  distribution-  with  the  theoretical  steady  state 
distribution  at  various  times  from  reservoir  filling  (early 
1973)  to  the  end  of  1978*  This  should  indicate  whether  or 
not  the  pore  pressure  response  due  to  filling  and  transient 
flow  is  converging  on  the  expected  steady  state  conditions# 

A  number  of  different  methods  of  comparison  were 
adopted,  the  purpose  being  to  select  the  best  method  to 
visualize  the  pore  pressure  behaviour# 


5 .a 2.  g.ppjpftrlagfl  a* 


Sections 


£E24axa£S1 


Plots  of  piezometer  locations  and  of  piezometric  level 
versus  time  are  shown  on  Figs#5#l  to  5#4#  These  figures 
include  the  theoretical  steady  state  response  of  the 
piezometers#  During  transient  seepage  it  would  be  expected 
that  the  actual  readings  would  converge  on  the  steady  state 


readings# 

Some  trends  are  Immediately  apparent#  PE24  is  reading 
much  below,  and  does  not  tend  to  converge  on,  the  steady 
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FIGURE  5.1 


Piezometric  Level  ( PE24)  versus  Time 


TIME  —years 
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FIGURE  5.2 


Piezometric  Level  (PH1.PE35)  versus  Time 


TIME  —  years 


Reservoir  level 


FIGURE  5.3 


Piezometric  Level  (PH  2  ,  PH  6,  PH  7)  versus  Time 


71 


o 


I  I 


FIGURE  5.4 


Piezometric  Level  (PHIO,  PH  11 )  versus  Time 


Tl  ME  —years 


* 
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state  condition  throughout  the  5  years  shown*  PH2  is  reading 
very  high  at  all  times  and  does  not  converge  on  the  lower 
steady  state  values*  However,  PHI  appears  to  be  reading  very 
close  to  the  steady  state  conditions  from  early  1976 
onwards*  The  actual  readings  are  a  little  high  (approx*  5m) 
but  this  is  only  an  error  of  6%* 

An  interesting  observation  is  that  if  the  steady  state 
line  is  shifted  vertically  it  tends  to  match  the  actual 
conditions  in  a  similar  manner  to  PHI*  PHI  suggests  that 
steady  state  conditions  have  been  reached  by  mid-1976  and 
that  prior  to  this  time  there  must  be  transient  seepage*  The 
time  lag  should  be  noted  as  the  reservoir  is  raised  to  a 
peak,  after  the  reservoir  level  drops  again  the  piezometer 
continues  to  rise*  This  behaviour  may  be  explained  as 
transient  seepage*  Even  though  the  reservoir  level  is 
dropping  the  steady  state  value  for  a  particular  piezometer 
may  still  be  higher  than  the  actual  value,  so  the  piezometer 
will  continue  to  rise* 

The  large  differences  between  calculated  and  observed 
values  for  PE24  and  PH2  cannot  be  explained  by  anisotropy, 
n on— homo gene ity  nor  any  other  parameter  already 
investigated*  Figure  5*5  clearly  shows  the  unusual  pore 
pressure  distribution  at  various  time  intervals* 

5 .*2x2  FBfcj.txEBZ 

Plots  of  piezometric  level  versus  time  are  shown  in 
Fig«5*3*  The  steady  state  case  is  superimposed  on  these 
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®  PE24  PHI  PH2  ® 


FIGURE  5,5  Pore  Pressure  Distribution  (22  +  50)  across  Elev.  590-7  m 

with  Time 
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figures*  It  has  already  been  noted  that  either  PR6  or  PH7 
must  be  incorrect*  A  comparison  with  the  theoretical  steady 
state  core  indicates  that  PH6  is  in  fact  very  close  to 
steady  state  from  1977  onwards*  wherees  PH7  is  actually 
reading  above  steady  state*  Once  again  if  the  steady  state 
line  is  displaced  vertically  the  behaviour  is  very  similar 
in  nature  to  PH6*  PH6  again  illustrates  the  time  lag 
associated  with  transient  seepage*  As  for  PE24  and  PH2  it  is 
impossible  to  explain  why  PH7  is  reading  too  high  from  the 
parameters  previously  studied* 

3  and 

Plots  of  piezoraetric  level  versus  time  ere  shown  on 
Figs*5*2  and  5*4*  The  steady  state  case  is  superimposed  on 
these  figures*  It  is  apparent  that  PE35  is  reading  much 
lower  (47  m)  than  the  predicted  steady  state  values*  even 
though  this  unit  has  been  recalibrated  upwards  to  account 
for  the  zero  shift  discussed  earlier*  Both  PH10  and  PHll 
appear  to  be  reading  too  high  by  about  6  m  and  9  m, 
respectively*  PHlO  and  PHll  are  both  perhaps  close  enough  to 
steady  state  conditions  that  the  differences  noted  could  be 
explained  by  a  number  of  minor  effects  such  as  reading  pore 
air  pressure  rather  than  pore  water  pressure  as  discussed 
earlier*  However»  there  is  no  explanation  fo  the  behavior  of 
PE35*  Once  again  the  piezometers  on  the  downstream  side  of 
the  core  are  reading  high  whereas  the  piezometers  on  the 


upstream  side  are  very  low* 
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Drawings  were  not  prepared  for  the  remaining 
piezometers*  PH12  behaved  in  a  manner  very  similar  to  PHll* 
The  other  piezometers  such  as  PE38,  PP36*  PE44,  PP44,  PP45, 
PP46,PP47  were  not  compared  because  in  the  case  of 
PE38***PP36  and  PE44* • *PP44  the  readings  have  already  been 
noted  as  unreliable  ( Nussbaum  1S76).  I  he  complicating 
factors  for  the  piezometers  at  high  elevations  such  as 
PP45**«PP47  render  their  response  difficult  to  analyze* 

gsoslygiops 

In  an  attempt  to  define  a  common  trend  for  all  the 
piezometers  contour  plots  of  pore  pressure  were  made  on  a 
composite  section  for  various  reservoir  levels*  They 
indicate  a  definite  pore  pressure  bulge  on  the  downstream 
side  of  the  core*  Although  it  may  not  be  as  conclusive  due 
to  the  lack  of  functioning  or  reliable  piezometers  on  the 
upstream  side,  both  PE24  and  PE35  are  considerably  below  the 
theoretical  steady  state  value* 

This  unusual  pore  pressure  distribution  does  not  show 
any  convergence  on  the  theoretical  steady  state  distribution 
over  the  5  years  observed  to  late  1978*  However,  some 
individual  piezometers  sdch  as  PHI  and  PH6  appear  to  show 
that  steady  state  has  indeed  been  reached  as  early  as  1976 
and  1977*  The  other  piezometers  with  a  vertical  "correction" 
also  show  this  trend*  The  obvious  disparity  from  the  steady 
state  values  cannot  be  explained  by  any  of  the  parameters 
previously  studied  in  the  parametric  evaluation* 


I  i 
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An  alternative  theory  is  obviously  required  to  explain 
the  observed  pore  pressure  distribution*  This  is 
investigated  further  in  Chapter  6* 


CHAPTER  VI 


INVESTIGATION  OF  PORE  PRESSURE  ANOMALIES 

In.tggrfagtioa 

Some  observed  pore  pressures  in  the  core  are 
significantly  different  to  the  predicted  steady  state 
values*  The  intent  of  this  chapter  is  to  Investigate  these 
pore  pressure  anomalies  in  detail* 

The  pore  pressures*  upon  reservoir  impounding*  are 
altered  as  a  result  of  two  processes*  these  are: 

1*  Total  stress  Changes* 

2*  Transient  seepage  (which  includes  three 
sub-processes )* 

a*  Consolidation  and  swelling  (void  ratio 
changes  ) 

b*  Fluid  density  changes 
c*  Changing  degree  of  saturation 
It  is  very  difficult  to  distinguish  between  the 
processes  because  they  may  be  occurring 
simultaneously  within  the  core*  However  by  studying 
the  contribution  of  each  process  separately  it  is 
possible  to  evaluate  the  relative  importance  of 
each*  Figure  6*1  shows  schematically  the  processes 
occurring  in  a  simple  section* 


77 


78 


0  <  t  <c-o  t  =  c*o 


End-of-  construction 


Transient  state  Steady  state 


Note 

ut_  Increase  in  pore  pressure  due  to  reservoir  load 
Up -Increase  in  pore  pressure  due  to  advance  of  wetting  front 
uc_  Decrease  in  pore  pressure  due  to  consolidation 


FIGURE  6.1  Schematic  Diagram  of  Processes  Occurring  in  core 
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£-9.2  for?  Pressures  E&sjulllng  Is aa  Isiftl  Sirssa  £b&qggaj. 

The  changes  in  total  stress  which  are  of  concern  in 
this  section  are  those  resulting  from  reservoir  loading*  As 
the  reservoir  is  raised  and  lowered  there  are  increments  of 
total  stress  applied  and  removed  from  the  core*  There  will 
be  a  reaction  of  the  pore  pressures  to  these  increments  of 
total  stress*  To  determine  these  reactions  it  is  necessary 
to  calculate  the  magnitude  of  the  stress  changes  at 
representative  locations  within  the  core*  Then  with  these 
values  and  appropriate  pore  pressure  parameters  the  pore 
pressure  changes  within  the  core  for  various  reservoir 
increments  may  be  calculated* 


£a2»1  Mac! el  lo.s 


F  l.Vl  i  h  g 


The  proposed  model  assumes  an  impermeable  upstream  core 
face*  This  means  that  the  reservoir  load  is  applied  directly 
to  the  upstream  face  of  the  core,  and  there  is  no 
interference  of  the  pore  pressure  distribution  due  to 
transient  seepage* 

The  reservoir  filling  does  not  impose  a  total  stress 
increment  on  the  core  face  equivalent  to  the  hydrostatic 
pressure  increment*  The  reason  for  this  is  that  the  upstream 
shell  is  only  slightly  saturated  initially  and  so  the 
reservoir  impounding  will  cause  complete  saturation  of  the 
shell*  The  horizontal  effective  stress  increment  is  related 
to  the  vertical  effective  stress  increment  by  means  of  the 


coefficient  of  earth  pressure,  K*  Hence  a  calculation  of  the 
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change  In  vertical  total  stress  (due  to  the  change  in  degree 
of  saturation)?  coupled  with  the  hydrostatic  pore  pressure 
and  K  will  yield  a  horizontal  total  stress  increment* 

As  shown  in  Figure  6*2  the  applied  horizontal  load  is 
approximately  60%  of  the  hydrostatic  load  for  this  section* 
The  order  of  magnitude  of  this  result  is  substantiated  by 
the  field  results  reported  for  the  Gepatsch  dam  in  Austria 
( Schober  1967)*  The  increment  in  vertical  stress  is 
calculated  to  be  quite  small?  which  concurs  with  Nobari  and 
Duncan  (1972)* 

The  incremental  load  is  then  distributed  across  the 
section  in  proportion  to  the  assumed  shear  stress 
distribution*  So  that  on  the  downstream  side  of  the  core  the 
increment  of  horizontal  stress  will  be  smaller  than  that  on 
the  upstream  side*  This  too  has  been  observed  on  the 
Gepatsch  Dam*  By  assuming  a  shear  stress  distribution  it  was 
possible  to  calculate  the  horizontal  and  vertical  stress 
Increments  at  some  points  across  the  core*  Using  these 
values  and  an  assumed  original  stress  distribution  obtained 
from  Simmons  (1974)  and  Law  (1975)  it  was  possible  to  draw 
the  Mohr  circles  and  stress  paths  followed  during  reservoir 
loading*  This  in  turn  provides  the  increments  in  principal 
stresses  A<r,  and  ^^3  which  may  be  used  to  determine  the  pore 
pressure  changes?  using  Skempton’s  (1954)  equation*  The  Mohr 
circles  and  stress  paths  are  shown  on  Fig*6*3*  and  Flg*6*4* 
An  initial  assumption  was  that  the  major  principal  stress 


was  vertical* 


, 
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SHEAR  STRESS  DISTRIBUTION 


FIGURE  6.2  Reservoir  Loading  Model  and  Results 
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p _ MPa 

ELEMENT  C 


p _ MPa 

ELEMENT  D 

NOTES:  —  T.S  P  —  Total  stress  path 

—  E.S.  P  _  Effective  stress  path 

—  The  calculated  E.S.P  assumes  zero  excess 
construction  pore  water  pressure 


FIGURE  6.4  Stress  Paths 


' 
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The  principal  stress  increments  calculated  for  this 
section  at  elevation  590*7  m  are  of  the  same  order  of 
magnitude  as  those  calculated  for  the  Oroville  Dam  by  Nobari 
and  Duncan  ( 1972)  and  as  those  observed  in  the  Gepatsch  dam 
by  Schober  (1970)*  So  it  is  felt  that  this  model  provides  a 
reasonable  indication  of  the  pricipal  stress  incrementst  and 
rotation  of  the  principal  stresses*  A  much  more  refined 
analysis  could  be  made  in  a  similar  manner  to  that  of  Nobari 
and  Duncan  (1972),  taking  into  account  the  change  in 
stress— strain  characteristics  of  the  upstream  shell  on 
wetting  and  so  on*  But  the  pore  pressure  increments  are 
dependent  on  the  pore  pressure  parameters,  and  without  more 
accurate  pore  pressure  parameters,  the  refined  stress 
increments  would  be  largely  wasted  effort* 

6*2*2  Poig  Prgssara  Ea j&me.terg 

The  total  stress  increments  have  been  calculated  due  to 
the  reservoir  loading*  In  order  to  predict  pore  pressure 
changes  due  to  these  stress  increments  it  is  necessary  to 
use  a  relationship  of  the  form  suggested  by  Bishop  ( 1954), 
Skeropton  (  1954),  Henkel  (I960  )* 

Bishop  (1954)  suggested  the  use  of  the  parameter  B, 

where 


6.  1 
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Af  B  are  Skempton' s  pore  pressure  parameters* 

AOj  is  the  minor  total  principal  stress  increment* 
is  the  major  total  principal  stress  increment* 
Then  the  basic  pore  pressure  equation  becomes: 


6.2 


•  •  • 


AU  is  the  pore  pressure  increment* 


A<r,  is  the  major  total  principal  stress  increment* 

B  is  the  pore  pressure  parameter* 

The  advantage  of  using  this  parameter  B  is  that  the 
pore  pressure  can  be  expressed  as  a  direct  function  of  the 
major  principal  stress  increment*  Unfortunately,  B  is  a 
function  of  the  total  principal  stress  increment  ratio,  in 
cases  where  this  ratio  changes  significantly  the  B  concept 
holds  no  advantage  over  Skempton's  (1954)  equation* 


The  equation  suggested  by  Skempton  takes  the  form: 


Aa  =  +  (\CAer;- A^)J 


•  •  • 


6*3 
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The  pore  pressure  parameters  are  determined 
exper imen ta lly , ho we ver  A  is  stress  dependent  and  B  is  a 
function  of  the  degree  of  saturation  of  the  soil*  This 
equation  makes  no  allowance  for  the  intermediate  principal 
stress  because  it  has  been  derived  for  the  special  case  of 
the  triaxial  test  where  • 

In  I960*  Henkel  proposed  a  similar  but  more  general 
expression*  with  different  experimentally  derived  pore 
pressure  parameters.  Henkel's  equation  takes  the  form: 


6.  4 


where  o4.  is  the  pore  pressure  parameter  related  to  the 


octahedral  normal  stress  Increment  and 


is  the  parameter 


related  to  the  incremental  octahedral  shear  stress. 

This  equation  is  more  suitable  than  Skempton's  (1954) 
equation  in  that  it  better  defines  the  soil  behaviour. 
Particularly  for  the  complicated  stress  paths  followed  in  a 
dam  core.  To  use  this  equation  It  is  necessary  to  make  an 
assumption  about  •  This  may  be  done  by  assuming  plane 

strain.  A  valde  for  is  required.  ^  may  be  evaluated 

from  A  obtained  from  triaxial  testing.  The  major  advantage 
of  this  approach  is  that  ^  is  not  stress  path  dependant. 
The  O L  value  at  high  stress  levels  may  not  be  unity.  The 


value  of  Poisson's  ratio  used  was  0.35  (Law  1S76  )•  The  pore 
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pressure  Increments  calculated  by  this  method  are  compared 
to  the  values  calculated  using  Skempton* s  equation#  A 
problem  arises  in  the  use  of  Skempton's  equation  for  loading 
cases  other  than  the  standard  triaxial  compression  test* 
Skempton  (1961T  discussed  the  case  of  extension  loading  and 
introduced  an  absolute  value  of  total  principal  stress 
difference;  however,  an  ambiguity  still  exists  in  the 
definition  of  principal  stress  increments,  especially  in  the 
particular  case  where  there  is  rotation  of  the  principal 
st  resse  s • 

Law  and  Holtz  (1978)  suggest  an  approach  to  deal  with 
this  ambiguity  and  their  approach  has  been  adopted  for  this 
analysis*  In  effect,  they  redefine  and  as  the 

algebraic  maximum  principal  stress  increment  and  minimum 
principal  stress  Increment  respectively,  and  the  Increments 
have  no  specific  relationship  with  the  principal  stresses* 
This  Involves  a  slight  redefinition  of  the  A  parameter 
depending  on  the  specific  triaxial  test  under  consideration* 
They  note  that  the  A  value  for  an  extension  test  is 
approximately  double  that  for  a  compression  test,  if  the 
material  is  behaving  approximately  elastically*  They  suggest 
t  hat  this  system  removes  the  amblguiti es  which  deve lop  when 
the  principal  stresses  rotate*  However,  since  they  give  no 
specific  example  of  this  phenomenon  their  interpretation  may 
need  further  definition* 

The  rotation  of  the  principal  stresses  may  be 
equivalent  to  an  increase  in  A  as  reported  by  Brotns  and 
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Casbarian  ( 1965);  howevept  without  further  evidence  and  test 
results  for  a  compacted  till  this  effect  is  ignored  as 
insignificant  due  to  the  very  low  A  values  observed  during 
standard  testing* 

It  is  necessary  to  ascertain  a  value  for  A*  Insley  and 
Killis  (1965)  published  the  results  from  an  extensive  series 
of  high  pressure  triaxial  tests  carried  out  on  the  Mica 
till*  These  were  undrained  axial  compression  tests*  They 
show  a  plot  of  the  parameter  versus  *  where  A  equals 

AxB  (Bishop  and  Henkel*-  1962)*  This  is  not  what  was 
required,  since  the  core  will  not  even  be  close  to  failure 
during  reservoir  loading*  So  some  A  values  were 

backcalcula t ed  from  the  published  data  at  specific  confining 
pressures*  The  results  are  shown  on  Fig*  6*5*  They  indicate 
low  but  positive  values  of  A  at  the  strains  estimated  in  the 
core*  It  would  be  expected  that  the  A  value  should  be  close 
to  0*5  if  the  soil  is  behaving  elastically  in  plane  strain* 
The  low  observed  values  are  probably  a  result  of  incomplete 
saturation  in  the  samples  tested*  Some  results  from  a 
compacted  moraine  material  reported  by  Bishop  and  Henkel 
(1962)  indicate  a  negative  Ap  of  the  order  -0*1*  In  other 
words  the  A^  values  reported  by  Insley  and  Hillis  (1965)  for 
the  Mica  till  are  higher  than  expected  and  they  attribute 
this  behaviour  to  the  rather  higher  than  normal  confining 
pressures • 

The  A.  values  to  be  used  in  this  enalysis  vary  between 


0*04  and  0*08  for  axial  compression  loading  and  double  these 


' 
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FIGURE  6.5  "A 


Values  for  Mica  Till 
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values  for  lateral  compression  values*  depending  on  stress 
levels  at  the  element  being  analyzed* 

The  B  parameter  is  a  function  of  degree  of  saturation* 
The  degree  of  saturation  of  the  core  at  the 

end-of-cons truct 1  on  has  been  estimated  (Section  6*2*3)  to  be 
approximately  93%  at  elevation  590  m*  This  results  In  a  B 
value  of  approximately  0*7  ( Lambe  and  Whitman,  1969  — 

Boulder  clay)*  This  value  may  be  slightly  in  error,  but  it 
will  not  introduce  a  significant  error*  As  a  simplification 
it  was  assumed  that  the  B  value  will  increase  to  0*9  during 
reservoir  filling*  This  value  may  be  a  little  high  (Lee  et 
el*  1969)  however  the  results  are  not  very  sensitive  to  this 
var iat ion* 

6.2.3  Degree  a&tug&tjJUr  cfir.£ 

The  degree  of  saturation  of  the  core  till  at  the  time 
of  placement  can  be  calculated  knowing  the  specific  gravity, 
dry  density  and  bulk  density*  In  this  manner  an  average 
degree  of  saturation  was  evaluated  for  both  the  lower  (below 
elevation  610  m)  and  upper  zone  of  the  core*  Some  typical 
properties  are  tabulated  below: 


Property 

Uni  t 

1969  Fill 

1970-1971 

tfd 

Mg/  m3 

2.  10 

2.27 

w*  c  • 

% 

8.8 

6.3 

Yt 

Mg/  m3 

2.29 

2.41 

Gs 

2.8 

2.8 

Sr 

% 

74 

74.5 

' 
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n  0.250  0.190 

The  specific  gravity  of  the  Mica  till  was  reported  as 
2.80  (Law  1975).  This  value  appeared  high  so  some  specific 
gravity  tests  were  carried  out  which  confirmed  the  high 
values.  The  degree  of  saturation  will  increase  with  the 
addition  of  fill.  With  some  simplifying  assumptions  it  is 
possible  to  backcal cul at e  the  changing  degree  of  saturation 
as  construction  proceeds  ,  using  the  measured  increases  in 
pore  pressure.  The  method  is  based  on  that  of  Hilf  (1948)  or 
B ishop  (  1957  )• 

Calculations  for  PE24  and  PH6  indicate  a  final  degree 
of  saturation  in  the  lower  core,  before  reservoir  filling 
comences,  of  94  %• 

These  values  show  that  the  pore  pressure  parameter  B 
will  certainly  not  be  unity.  Larabe  and  Whitman  (1969) 
indicate  a  value  for  B  at  93  %  saturation,  for  a  boulder 
clay.  This  value  is  considered  to  be  reasonably 
representative  for  the  Mica  till. 

The  increase  in  pore  pressure  is  caused  by  a  volume 
change,  which  is  initiated  by  the  increase  in  total  stresses 
on  the  unsaturated  fill.  The  high  compressibility  of  the  air 
governs  the  volume  change  and  the  diffusion  of  the  air  into 
the  pore  water  due  to  the  increase  in  pore  pressure.  The 
laws  which  govern  this  behaviour  are  Boyles  Law  and  Henry's 
Law  of  Solubility.  The  volume  changes  calculated  due  to  the 
changes  in  effective  stress  should  be  compatible  with  the 
volume  change  calculated  by  the  method  outlined  above. 
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The  principal  stress  increments  have  been  calculated 
for  various  stages  of  reservoir  filling*  It  was  decided  to 
use  Skemptons  (1954)  pore  pressure  equation  to  evaluate  the 
pore  pressure  increments*  Values  for  A  and  B  have  been 
assumed*  Therefore  using  the  equation: 


6.5 


the  pore  pressures  may  be  calculated*  As  discussed  in 
Section  6*2*2  the  approach  suggested  by  Law  and  Holtz  (1978) 
is  adopted  here,  for  determining  the  maximum  and  minimum 
principal  stress  increments* 

The  results  of  this  analysis  are  shown  on  Figure  6*6. 
The  A  value  has  been  altered  with  each  reservoir  stage.  Some 
allowance  has  therefore  been  made  for  stress  level  effects 
on  A,  that  is,  the  reduction  of  stresses  in  the  central  core 
due  to  arching*  The  alteration  in  the  A  value  depending  on 
the  stress  path  has  been  mentioned  previously  and  this 
factor  is  also  taken  into  account*  During  the  early  loading 
the  stress  path  is  similar  to  that  of  an  ordinary  triaxial 
extension  (lateral  compression)  test  and  later  this  changes 
to  that  of  a  triaxial  compression  (axial  compression)  test* 
Which  is  why  the  A  value  has  been  reduced  for  the  second  and 
third  loading  stages*  The  final  factor  is  that  the  A  value 
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is  dependent  on  the  mobilized  strength  cf  the  element'  and 
this  factor  has  also  been  estimated*  A  combination  of  these 
factors  result  in  the  selected  values  for  A  between  0*05  and 
0.  10. 

The  B  value  has  simply  been  assumed  to  Increase  with 
the  increasing  loading*  This  may  be  an  over-simplification 
and  the  final  value  of  0*9  is  probably  too  high*  however  it 
has  little  effect  on  the  results* 

2.8.5-  gyi.labj.lj  ly  Scooted  ErQQ.S^.iAlLSs 

It  is  possible  that  doubt  may  be  expressed  over  the 
approach  taken  by  Law  and  Holtz  (1978)*  The  deviation  from 
the  standard  approach  only  occurs  during  the  first  loading 
stage*  The  pore  pressure  increments  calculated  from  this 
stage  were  compared  with  the  measured  pore  pressure 
increments  calculated  from  Henkel's  (1960)  equation*  These 
are  tabulated  below: 

Element  Res*  Incr*  Henkel(  1960)  Sk empt on(  1 9 54  ) 

Au  —  m*  Au  — m* 

B  1  8*5  9*8 

2  12.0  14.3 

3  9.3  8.8 

Total  29.8  32.9 

C  1  6.6  7*1 

2  4.4  3.2 

3  3.8  1.1 

14.8  11.4 


Total 
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D 


1 


4.4 


6.6 


2 


3.9 


6.6 


3 


0.2 


-2.5 


Total 


8.1 


7.7 


The  comparison  is  excellent  considering  the  assumptions  used 
in  the  analysis.  This  appears  to  confirm  that  Skempton's 
( 1954)  equation  may  be  used  provided  care  is  taken  in  the 
definition  of  the  principal  stress  increment.  Henkel's 
(1960)  equation  is  more  suitable  for  complicated  stress  path 
problems . 

The  total  stress  /  pore  pressure  increments  are  shown 
on  Figure  6.6.  The  conparispn  with  the  actual  pore  pressures 
is  excellent  at  the  early  stage  in  filling  (August  1973). 
There  is  little  effect  of  transient  seepage  so  the 
comparison  is  straightforward.  If  the  normal  procedure  had 
been  adopted  for  the  principal  stress  increments  in 
Skempton's  (1954)  equation  then  the  calculated  pore  pressure 
Increments  would  have  been  much  higher  than  the  observed 
values • 

Since  the  results  for  this  early  period  appear  to  be 
confirmed  by  the  observed-  values*  it  is  felt  that  the  later 
values  of  the  pore  pressure  calculated  by  the  same  method 
are  also  realistic.  It  is  difficult  to  check  these  due  to 
the  unknown  contribution  of  transient  seepage. 

Figure  6.6  shows  the  resulting  pore  pressure 
distribution  for  the  final  reservoir  level*  assuming  no 


further  consolidation  on  the  downstream  side  of  the  core 
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which  is  a  conservative  assumption*  The  steady  state 
condition  and  the  actual  readings  are  shown*  It  would  appear 
then  that  the  pore  pressures  induced  due  to  reservoir 
filling  are  not  sufficient  to  explain  the  observed 
d is  t  ributl on • 

One  interesting  result  of  the  analysis  is  that  the 
element  close  to  the  downstream  face  actually  exhibits  a 
slight  negative  pore  pressure  increment  during  the  final 
stage  of  reservoir  filling*  This  phenomenon  is  also  apparent 
in  Chang’s  (1976)  theoretical  study  of  the  New  Melones  Dam 
(California)*  Even  though  he  uses  a  completely  different 
approach  for  his  analysis  ( Caro— clay  model),  he  shows  in  some 
zones  of  the  core  a  reduction  of  pore  pressure  on  reservoir 
loading* 


6a  3  Izagflull-C  Era  slurs 

Hydraulic  fracturing  within  embankment  dams  is  a 
process  which  has  only  recently  been  recognised  by 
Geotechnical  Engineers*  Kjaernsli  and  Torblaa  (1968) 
suggested  this  mechanism  as  a  possible  cause  of  leakage 
through  the  Hyttejuvet  Dam*  Sherard  (1973)  provides  a  good 
review  of  some  notable  case  histories  which  exhibited  this 
phenomenon*  Since  this  time  the  process  has  been  observed  in 
other  dams,  and  has  even  been  purposely  initiated  from 
boreholes  within  a  dam  core  (Penman  1975)*  A  recent  paper  by 
Kulhawy  and  Gurtowski  (1976)  studied  the  parameters  likely 
to  affect  hydraulic  fracture  potential  through  a  series  of 
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finite  element  analyses*  The  volume  of  discussion  arising 
from  this  paper  indicates  the  interest  presently  held  in 
this  area*  Hydraulic  fracturing  is  now  a  widely  recognised 
phenomenon  and  embankment  designers  must  consider  this 
process  in  their  design*  To  date  there  is  an  incomplete 
under standi ng  of  the  mechanisms  Involved  end  disagreement  as 
to  which  parameters  increase  the  likelihood  of  hydraulic 
fracture  (Charles  1976)* 

This  section  analyzes  the  likelihood  of  hydraulic 
fracture  occurring  in  the  Mica  Dam*  Use  is  made  of  the  total 
stress  increments  calculated  earlier  in  Section  6*3* 

Hydraulic  fracturing  has  allegedly  occurred  in  the  Mica 
core  (Simmons  1974,  Skermer  1975,  Nuesbaum  and  Stevenson 
1976)  during  inflow  tests  on  the  vertical  movement  gauge 
MV15*  The  test  was  carried  out  when  the  reservoir  was  at 
elevation  676  is*,  the  maximum  level  to  which  the  water  could 
be  raised  in  the  gauge  was  666  m*  Both  Simmons  and  Skermer 
show  that  the  hydraulic  fracture  could  have  Initiated  below' 
elevation  590  m*  However  neither  considered  the  increase  in 
the  minor  principal  stress  due  to  the  reservoir  filling  to 
elevation  676  in*  The  total  stress  change  analysis  carried 
out  in  this  thesis  indicated  a  change  in  minor  principal 
total  stress  of  about  408  kPa  (3*S  tsf  )  which  would  mean 
that  the  hydraulic  fracture  could  not  occur  above  elevation 
533  m  which  is  about  5  m  above  the  bedrock*  This  assumes 
that  the  results  of  Simmons  ( 1974)  3— Dimensional  finite 


element  analysis  are  correct*  The  condition  of  raising  the 
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water  level  In  a  borehole  or  vertical  movement  gauge  however 
is  not  the  same  as  raising  the  reservoir*  Raising  the 
reservoir  causes  total  stress  increments  within  the  core 
which  should  not  be  ignored*  If  the  dam  is  arched  between 
the  abutments  the  stresses  may  be  redistributed*  The  effect 
of  wetting  the  upstream  shell  may  have  a  further  influence 
on  the  stresses  so  it  cannot  be  evaluated  with  certainty 
whether  hydraulic  fracturing  is  likely  or  reservoir  filling 
from  the  results  of  a  simple  analysis  or  from  a  simple 
inflow  test* 

Hydraulic  fracture  potential  must  be  evaluated  on  the 
upstream  face  of  the  core*  The  stress  state  at  this  location 
is  difficult  to  establish  accurately  due  to  the  large  stress 
gradients  across  the  interface*  If  it  is  assumed  that  the 
core  arching  is  continuous  to  the  edge  of  the  core»  then 
using  Simmons  (1974)  stress  results  for  elevation  590  m* , 
the  minor  principal  stress  is  850  kPa  (8  tsf)»  whereas  the 
hydrostatic  pressure  at  this  level  is  1800  kPa  (17  tsf  )• 
However  this  does  not  include  the  increase  due  to  the 
reservoir  loading*  This  is  estimated  to  be  850  kPa  (8  tsf)* 
therefore  the  minor  principal  stress  is  actually  about  1700 
kPa*  This  is  very  similar  in  magnitude  to  the  reservoir 
pressure, and  if  the  values  are  accurate  then  hydraulic 
fracture  would  probably  not  occur  as  it  seems  to  initiate  at 
some  value  between  and  (Penman  1972),  this  may  be 

due  to  some  structural  behaviour  or  simply  because  the  soil 
has  some  tensile  strength*  An  added  consideration  is  that 
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the  direction  of  C~j  is  not  conducive  to  fracture 
development,  and  any  fractures  which  develop  are  likely  to 
form  perpendicular  to  (vertical  cracks)  or  in  some 

direction  perpendicular  to  an  average  of  0^  and 
(horizontal  or  sub-horizontal  cracks)*  The  major  principal 
total  stress  calculated  (including  reservoir  load)  is  257S 
kPa  compared  to  a  reservoir  pressure  of  1800  kPa<  Therefore 
horizontal  or  sub— horl zon tal  cracks  are  not  likely  to  occur 
in  this  case* 

The  results  do  indicate  that  if  significant  arching  (60 
%  or  greater)  occurs  in  the  core,  and  this  arching  is 
present  to  the  edge  of  the  core,  then  hydaulic  fracture  is 
indeed  possible  even  though  the  reservoir  load  does  increase 
the  principal  stresses  in  the  core* 

The  use  of  controlled  hydraulic  fracture  tests  in  dams 
may  prove  very  useful  in  the  future  to  evaluate  the 
likelihood  of  hydraulic  fracture  occurring  during  reservoir 
filling*  This  method  would  allow  the  designer  to  evaluate 
the  degree  of  arching  occurring  before  reservoir  filling  and 
hence  a  more  accurate  feeling  for  the  fracture  potential* 
However  much  more  research  work  is  required  in  this  area  to 
provide  insight  into  the  physical  processes  involved* 


Seepage 


Transient  seepage  is  the  second  major  process  occurring 
in  the  core  during  resevoir  impounding*  A  governing  equation 
is  required  which  includes  the  three  sub-processes  occurring 
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within  the  transient  state*  This  flow  equation  should 
provide  the  pressure  head  distribution  with  time  for 
transient  flow  through  an  unsaturated  porous  medium* 

Sgyer, nlpg  Eauaiioa 

The  equation  of  continuity  for  transient  unsaturated 
flow  Is  (Freeze  and  Cherry,  1979)  : 


6*6 


f  -  Mass  density 
v  =  velocity 
n  =  porosity 

s  =  degree  of  saturation 
Following  from  this  the  equation  of  flow  is  : 


■*x 


+ ft  0] .  0^) 


6*7 


Pressure  head 


k  =  Permeability 


C  =  Specific  Moisture  Capacity 
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This  equation  is  often  known  as  Richards  equation*  The 
solution  for  (x,y,z,t)  describes  the  pressure  head  field  at 
any  point  in  a  flow  field  at  any  time* 

The  specific  moisture  capacity  in  equation  6*7  is 
further  defined  as  : 


Sn  XT  ±  S  yn 
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6.8 


where  0  =  Volumetric  moisture  content 

The  term  is  the  term  relating  fluid  density  changes* 

J* 

The  term  S  i s  the  term  which  considers  void  ratio  changes, 

the  last  term  is  the  change  in  degree  of  saturation* 

yf 

Normally  the  fluid  density  and  void  ratio  change  terms 
are  ignored  for  unsaturated  soil  problems  (Freeze  and 
Cherry,  1979  )  because  they  -are  small  with  respect  to  the 


change  in  degree  of  saturation  term* 

This  is  usually  true  for  the 
consolidation  or  swelling  term  should  be  investigated 


5*?  V  term,  however  the 

r  *** 


further  to  determine  its  magnitude  for  each  case* 

If  the  soil  considered  is  near  full  saturation  then  the 


r elat ionsh i p 


<r '  =■ 


<3~~  —  L/< 


9  9  9 


6.  9 


102 


will  apply  where 

1 

O  =  effective  stress 
G~~  =  total  stress 

The  term  may  be  redefined  as 


-S 


s  >£  .  1°: 
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Now  where  ^v'  is  the  coefficient  of  volume 


change • 


22  =  Vw . 
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Therefore  for  a  soil  near  full  saturation  the  specific 
moisture  capacity  becomes  approximately 


c<y)  = 


6.  13 


When  this  is  replaced  in  eqn,  6,7  then  the  classical 
consolidation  equation  is  recovered,  assuming  that  K»  the 
permeability,  is  unique  within  the  pressure  range  under 
con3iderat ion. 

For  the  case  where  there  is  incomplete  saturation  then 
equation  6,7  is  not  sufficient  to  describe  the  behaviour.  In 
order  to  compare  the  magnitude  of  the  terms  S and 


upper  bound  value  for  S 


Using  this  assumption  the  specific  moisture  capacity 
then  becomes 


6,14 


For  practical  purposes  the  situation  where  equation 


6,14  may  be  of  concern  is  where  is  increasing  and  the 


first  term  then  relates  to  swelling. 


Hence 
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CMs)  -  _S  Yw  -h  n  *5 


6.15 


where  /v>^  coeff.  of  volume  change  w.r.t.  swelling. 


A  calculation  done  on  the  Mica  dam  core  material  indicates 


an  order  of  magnitude  difference  between  the  first  and 
second  terras  of  equation  6.15.  Hence  i t  is  sufficiently 
accurate  to  ignore  the  swelling  term  in  the  governing 
equation  and  an  analysis  considering  a  rigid  skeleton  should 
yield  sufficiently  accurate  results  in  the  zone  of  the  core 
where  the  pressure  is  increasing. 

An  error  may  be  involved  in  the  downstream  zone  of  the 
core  where  consolidation  is  continuing.  To  avoid  the 
complication  of  ascertaining  the  consolidation  /  swelling 
boundary  required  for  an  accurate  solution  of  equation  6.7, 
an  as sump ti on  was  mode  that  t  here  is  full  sa tura  ti on  across 
the  core  and  further  that  approximately  equals  and  the 

resulting  consolidation  equation  is  one  dimensional  nature. 


This  results  in  a  specific  moisture  capacity  of 


•  .  .  6.16 


Hence  equation  1  becomes 
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6.  17 


where  =  coefficient  of  consolidation. 

This  is  Terzaghl 1 s  one  dimensional  consolidation  equation. 
It  is  simply  solved  by  a  finite  difference  explicit  sc  heme, 
otherwise  known  as  the  marching  forward  technique. 


The  resevolr  impounding  causes  an  increase  in  pressure 
head  on  the  upstream  face  of  the  coret  which  is  transmitted 
throughthe  core  by  the  process  of  transient  seepage.  The 
assumption  of  a  rigid  soil  skeleton  is  common  for  transient 
flow  in  unsaturated  porous  media  (Freeze  and  Cherryf  1979). 
In  effect  what  this  means  is  that  both  of  the  terms  in  the 
specific  moisture  capacity  equation  6.8»  relating  to  void 
ratio  changes  and  fluid  density  changes  are  negligable  in 
comparison  to  "the  the  saturation  term.  In  other  words 


•  •  • 


6.18 


or  in  an  alternative  form 


•  •  • 


6.19 
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This  validity  of  this  assumption  was  checked  for  the  Mica 
till  using  an  hi ^  value  equivalent  to  one  fifth  the 
value.  This  resulted  in  an  order  of  magnitude  difference 
between  Syw  Mj  and  n  which  supports  the  assumption  of  a 

rigid  soil  skeleton. 

A  review  of  the  literature  provided  a  number  of 
solutions  to  the  transient  flow  problem*  with  the  assumption 
of  a  rigid  soil  skeleton.  These  were  by  Brahma  and 
Harr(  1962  )  ,  Desai  and  Sherman  (1971)*  Dvinof f(  1970  )  ,  Dvlnoff 
and  Harr(1971)r  and  Stephenson( 1978  )•  For  this  thesis  the 
approach  suggested  by  Dvi nof f ( 1 970 )  was  followed. 

Dvlnoff  uses  the  linearized  governing  equation: 


•  •  • 


6.20 


h  =  total  head  at  a  point 
H  =  mean  height  of  reservoir 
N  =  Specific  moisture  capacity  x  H 
k  =  horizontal  permeability 
k  =  vertical  permeability 
x  =  horizontal  coordinate  direction 
y  =  vertical  coordinate  direction 
t  =  time 


The  assumptions  inherent  in  this  solution  are: 
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—  flow  Is  two  dimensional 

—  inertia  forces  are  insignificant 

—  fluid  is  continuous 

—  air  escapes  freely  from  pore  spaces 

—  soil  is  rigid,  n,k  are  constant  with  time 
—  Darcys  law  is  valid 

—  DupuLt's  assumptions  are  valid 

In  order  to  use  the  charts  provided  by  Dvinoff  (1970), 
some  simplifying  assumptions  were  made  regarding  the 
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analyses  to  actual  case  histories,  and  so  have  not  indicated 
typical  values  for  this  'storage*  coefficient,  save  for 
their  model  tests  where  N  is  simply  the  total  porosity* 

It  is  assumed  for  this  study  that  N  may  be  evaluated 
from  the  degree  of  saturation  existing  at  the  end  of 
construction*  That  is,  assuming  a  rigid  soil  structure 
calculate  the  change  in  degree  of  saturation  required  to 
achieve  a  steady  state  pore  pressure  distribution,  this 
coupled  with  the  actual  porosity  will  give  the  required 
volume  of  inflow  per  unit  volume  (storage  coefficient)*  This 
value  will  not  be  the  same  under  drawdown  conditions  because 
of  capilliary  effects*  However  the  reservoir-time  curve 
assumed  ddes  not  include  any  drawdown*  The  storage 


* 
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Linear  approximation  to 
1  reservoir  level 
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FIGURE  6 .7  Reservoir  Filling  Curve  and  Approximations 
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coefficient  will  not  necessarily  be  constant  across  the  dam 
so  an  average  value  has  been  chosen  from  the  central  core 
zone  • 

The  approach  adopted  for  solution  of  the  composite 
reservoir— time  curve  is  similar  to  that  of  Dvinof f  (  1970)* 
That  is,  simplified  paths  are  assumed  as  upper  and  lower 
bounding  conditions,  each  of  which  may  be  solved  using  the 
charts  and  tables  provided*  Then  the  composite  solution  is 
approximated  between  these  limiting  solutions*  The  results 
are  shewn  on  Figure  6*8  for  the  time  intervals  of  1,  2  and  3 

years  after  reservoir  filling  commences* 

The  piezometers  were  plotted  on  these  figures  and  the 
calculated  piezometric  levels  were  then  plotted  versus  time* 
Superimposed  on  t  hese  plots  a  re  the  actual  piezometer 
readings  and  the  steady  state  values*  These  results  make  no 
allowance  for  pore  pressure  changes  due  to  the  reservoir 
loading*  These  effects  have  been  shown  to  be  relatively 
small  and  so  their  contribution  would  be  to  speed  up  the 
advance  of  the  phreatic  surface  by  a  small  ammount  (  the 
storage  coefficient  is  slightly  reduced)* 

6.4*3  jjj  as us si  on  qJL  r<?su.lia  £jl om  ir.anfele.al  sfeee&fce  sjia.iyg.lfe 

The  results  of  this  analysis  show  that  the  phreatic 
surface  proceeds  rapidly  through  the  core*  even  though  the 
permeability  is  very  low*  This  is  because  the  advance  is  not 
a  function  of  permeability  alone,  but  of  the  square  root  of 
permeability  divided  by  a  storage  coefficient  which  is 
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FIGURE  6. 8  Wetting  Fronts  for  Various  Reservoir  Paths 
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related  to  the  specific  moisture  capacity*  Errors  in  this 
function  are  reduced  by  virtue  of  the  square  root  term*  The 
value  for  the  specific  moisture  capacity  used  in  this 
analysis  was  0*0002  m— 1 • 

Evi dence  to  support  the  the ore  t ical  mo vemen  t  of  the 
wetting  surface  may  be  seen  in  the  observed  water  levels  in 
the  vertical  and  near— vert ica 1  movement  gauges*  These  values 
are  shown  in  Figure  6*9  and  compared  with  the  phreatic 
surface*  The  comparison  is  good  for  liV8  and  MV  1  5  »  however 
the  results  for  MV9  are  unusual*  The  water  level  begins  to 
raise  after  one  year  which  corresponds  well  with  the 
calculated  data*  However  the  value  at  the  end  of  2  years 
appears  rather  low*  This  could  be  due  to  inaccuracies  in  the 
solution  towards  the  base  of  the  coret  but  the  gauge  rises 
only  marginally  over  the  next  3  years*  This  behaviour  makes 
the  readings  of  MV9  suspect*  The  cause  could  be  some 
irregularity  in  construction  procedure  allowing  preferred 
seepage  path  to  the  downstream  shell.An  alternative  solution 
is  proposed  by  Freeze  (1971)  who  suggests  that  present 
steady  state  analyses  which  only  consider  flow  within  the 
saturated  zone  result  in  an  erroneous  phreatic  surface*  The 
complex  problem  of  assuming  flow  within  the  unsaturated  zone 
above  the  phreatic  surface  cannot  be  modelled  by  the 
simplified  transient  flow  analysis  used  here*  Hence  an  error 
in  the  solution  close  to  the  downstream  boundary  is 
possible* 


The  values  Indicated  by  MV's  4  and  14  show  a  more  rapid 
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FIGURE  6.9  Wetting  Front  Compared  with  Movement  Gauge  Readings 
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movement  of  the  phreatic  surface,  however  these  movement 
gauges  are  "located  in  much  shallower  sections  closer  to  the 
abutments  and  so  their  behaviour  may  not  necessarily  be  the 
same  as  for  the  deeper  section* 

More  interesting  are  the  plots  of  the  piezometric 
levels  versus  time,  Figs*6*10  and  6*11*  These  show  quite 
clearly  the  rapid  development  of  'approximate  '  steady  state 
conditions*  A  comparison  of  the  theoretical  values  with  the 
observed  values  may  be  made*  The  agreement  between  PHI,  PH6 
and  PE35  are  remarkably  good  considering  the  simplifications 
necessary  for  the  analysis*  The  earlier  observed  reactions 
of  these  piezometers  could  be  due  to  the  assumption  of  no 
pore  pressure  change  due  to  reservoir  filling*  The  general 
shape  of  the  reaction  and  the  convergence  cn  steady  state 
values  is  excellent*  However  large  inconsistencies  exist  in 
the  remaining  piezometers  PE24,  PH2 ,  PH7  and  PH10*  In  order 
to  evaluate  their  responses  it  is  necessary  to  consider 
their  respective  locations  (Figure  2*3  )• 

In  the  case  of  PE24  (upstream),  PHI  (central)  and  PH2 
(downstream),  PE24  reacts  at  the  correct  time  but  its 
response  is  much  too  subdued*  PHI  reacts  early  and  its 
response  is  correct*  PH2  reacts  much  earlier  than  the  theory 
suggests  and  its  reaction  is  above  the  theoretical  values 
for  steady  state*  Similar  behaviour  is  noted  for  PH6 
(central)  and  PH7  (downstream),  and  again  for  PE35 
(upstream)  and  PH10  (downstream)*  In  general  it  could  be 
stated  that  the  upstream  piezometers  react  at  the  correct 
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FIGURE  610  Advance  of  Wetting  Front  and  Calculated 

Piezometer  Readings 
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FIGURE  6.11  Predicted  Piezometric  Levels  during  Transient  Seepage 
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time  (approximately)  but  "their  response  is  low*  The  central 
piezometers  react  a  little  too  soon  and  their  response  is 
correct*  The  downstream  piezometers  react  much  too  soon  and 
their  response  is  too  high  (compared  to  steady  state 
theoretical  values)* 

In  general  it  is  considered  that  this  analysis  provides 
a  fairly  accurate  description  of  the  advance  of  the  phreatic 
surface  within  the  Mica  Dam  core*  The  advance  is  rapid  and 
is  substantiated  to  some  extent  by  the  pore  pressure 
readings  taken  in  the  vertical  ana  near— ver ti cal  movement 
gauge  s • 

The  assumption  that  the  consolidation  term  S ^5  in 
equation  6*8  is  negligible  compared  to  the  saturation  term 
may  not  be  valid  in  the  core  zone  where  the  pressure 

ty/ 

head  is  continually  reducing*  Particularly  if  the  degree  of 
saturation  is  high* 

In  order  to  investigate  the  rapidity  of  dissipation  of 
the  pore  pressures  in  the  consolidating  downstream  zone  the 
term  was  assumed  to  be  zero*  This  is  true  if  S=1  and  is 

not  an  unreasonable  assumption  for  this  material  at 
approximately  85%  saturation* 

To  simplify  the  calculations  the  upstream  ’swelling* 
zone  of  increasing  pressure  head  was  assumed  to  have  a 
coefficient  similar  in  magnitude  to  that  in  the 
consolidating  zone*  The  errors  Inherent  in  this  solution 
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will  occur  in  the  time  rate  of  advance  of  the  phreatic 
surface  in  the  upstream  zone*  However  the  indicated  rate  of 
consolidation  should  be  reasonably  accurate  for  the 
downstream  zone* 

The  coefficient  of  consolidation  was  assumed  to  be  14 
tn2/rao*  (Law  1P76)*  The  equation  was  solved  using  a  simple 
finite  di iference  (marching  forward)  explicit  1  Dimensional 
sc.heme.The  reservoir  Impounding  was  modelled  by  a  moving 
upstream  boundary  condition*  The  reservoir  level  was  held 
constant  at  elevation  734m  between  the  2nd  and  10th  years* 

The  solution  is  shown  on  Figure  6*12  •  The  results  show 
rapid  dissipation  close  to  the  downstream  face  of  the  core 
(PH2)»  with  only  minor  fluctuations  in  pore  pressure  as  the 
reservoir  rises*  Similarly  in  the  central  core  (  PHI  )  there 
is  little  decrease  in  pore  water  pressure  due  to 
consolidation*  It  can  clearly  be  seen  that  the  behaviour  of 
PH2  is  very  dissimilar  in  magnitude  to  the  theoretical 
behaviour.  During  the  early  stages  of  consolidation  however 
the  comparison  is  not  unreasonable  ( PB2  )•  It  is  not  expected 
that  the  upstream  side  of  the  core  ( PE24  )  will  match  well 
because  of  the  simplifications  introduced  into  the  governing 
equat ion* 

6*5  Cgnclusigaa 

This  chapter  reviews  the  physical  processes  occurring 
in  the  core  of  the  Mica  dam  prior  to  the  establishment  of 
steady  state  seepage*  The  processes  studied  separately  were: 
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FIGURE  6.12  Pressure  Head  Versus  Time 
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—  total  stress  changes  due  to  reservoir  filling 

—  transient  seepage 

The  most  significant  contribution  to  the  pore  pressu 
increments  was  found  to  occur  from  the  transient  flo 
advance  of  the  phreatic  surface  is  rapid,  approxinat 
4  years  to  develop  steady  stale*  Eecause  this  surfac 
nearly  vertical  the  response  of  the  piezometers  is  a 
rapid*  The  most  significant  contribution  to  the  tran 
process  is  from  the  changing  degree  of  saturation  fo 
unsaturated  soils  and  so  it  is  reasonable  to  essume 
soil  skeleton*  However  the  consolidating  portion  of 
is  more  affected  by  the  changing  void  ratio*  Whereas 
fluid  density  changes  are  not  normally  significant  a 
be  ignored* 

The  contribution  from  the  total  stress  changes 
undra in ed  re servoir  load i ng  was  calculated  to  be  sma 
not  negligible*  If  this  effect  was  to  be  superimpose 
development  of  the  phreatic  surface  it  would  reduce 
•storage1  coefficient  slightly  and  so  increase  the  r 
advance  of  the  phreatic  surface*  These  processes  are 
sufficient  to  explain  the  observed  pore  pressure 
distribution  in  this  embankment  core* 

A  comparison  of  the  observed  and  predicted  pore 
pressure  distributions  during  transient  seepage  or  s 
state  infers  some  irregularity  within  the  core*  A  po 
explanation  would  be  the  existence  of  hydraulic  frac 
the  core*  The  likelihood  of  fracturing  has  been  inve 
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and  the  results  indicate  that  hydraulic  fracture  is  unlikely 
in  the  upstream  zone  of  the  core,  but  would  be  possible  in 
the  central  and  downstream  core  zones  if  the  full  reservoir 
pressure  could  be  transmitted  in  some  manner  to  these  zones* 
However  because  the  unusual  pore  pressure  distribution  is 
similar  through  various  sections  and  levels  in  the  core  it 
could  be  argued  that  hydraulic  fractures,  if  they  did  occur, 
would  not  appear  as  uniformly,  nor  in  sufficient  quantity, 
to  create  this  distribution*  An  alternative  explanation 
could  be  that  some  material  non— homogenel t y  in  the  core 
construction,  yielding  much  more  impermeable  strips  parallel 
to  the  upstream  and  downstream  faces  of  the  core*  There  is 
no  justifiable  reason  to  assume  this  condition  would  be 
regularly  repeated  throughout  the  core,  and  there  is  no 
evidence  to  support  this  theory* 

None  of  the  theories  investigated  yield  results  even 
close  to  the  piezometric  values  observed*  Therefore 
consideration  must  be  given  to  the  fact  that  the  general 
pore  pressure  distribution  in  the  core  could  be  quite 
similar  to  the  theoretical  values  postulated,  but  for  some 
reason  there  are  some  local  anomalies  in  the  pore  pressure 
distribution  close  to  the  piezometers  and  so  these 
piezometric  readings  are  only  reflecting  local 
per t urba t ions* 

In  light  of  this  argument  the  reaction  of  piezometers 
PE38,  PP36  and  PE44  are  considered*  These  piezometers  are 
located  in  the  upper  core  and  are  noted  to  respond  almost 
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instantly  to  the  reservoir  fluctuation 
level  very  closely.  This  behaviour  has 
by  Nussbaum  (1976).  The  reason  for  thi 
is  because  the  instrument  lead  trenche 
directly  into  the  upstream  shell.  This 
although  the  remaining  trenches  do  not 
shell  they  may  have  a  local  influence 
distribution.  This  theory  is  investiga 
Chapter  7. 
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CHAPTER  7 


EFFECT  OF  PIEZOMETER  INSTALLATION  METHOD  ON  PORE  PRESSURE 

DISTRIBUTION 


In  recent  years  monitoring  o±  embankment  performance 
has  become  an  integral  component  of  the  design  procedure. 

The  reasons  for  this  are  twofold.  In  order  that  the  observed 
behaviour  can  be  compared  with  the  predicted  behaviour  and 
so  that  insight  into  the  physical  processes  occuring  during 
various  phases  in  the  life  of  the  embankment  may  be. 
devel oped . 

The  role  played  by  piezometers  is  extremely  important 
because  without  them  the  actual  pore  pressure  magnitude  and 
distribution  would  be  unknown*  excluding  the  possibility  of 
carrying  out  an  accurate  effective  stress  analysis. 

Three  basic  types  of  piezometer  are  in  common  use 
today,  each  having  distinct  advantages  over  the  others  in 
particular  situations.  However  no  matter  how  sophisticated 
the  instrument  the  readings  will  be  valueless  if  the 
installation  procedure  is  unsatisfactory.  The  one  common 
weakness  with  any  piezometric  system  is  the  link  between 
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measuring  unit  and  readout  unit*  Concern  has  been  expressed 

by  authors  that  these  lead  trenches  may  act  as  preferential 

seepage  paths  ( Sherard  et  al*  1963)*  He  states: 

" • • •  But  regardless  of  the  location  of  the  measuring 
point  f  the  tubing  must  be  carefully  backfilled  with 
Impervious  soil  to  prevent  seepage  along  the  line  of 
the  tubes  which  could  change  the  pore  water  pressure 
at  the  piezometer  tip*  " 

In  particular  Blight  (  1970)  noted  that  narrow  cores 
accentuate  this  danger,  and  considerable  attention  should  be 
given  to  design  and  installation  in  these  cases*  for  example 
he  quotes  the  Bridle  Drift  dam  (  S*  Africa)  and  the 
Man jiren ji  dam  (Rhodesia)* 

If  piezometers  are  installed  in  an  inadequate  manner 
the  observed  values  of  pore  pressure  may  be  misleading,  and 
in  certain  cases  lead  to  an  o  veres  tio.e  ti  on  of  the  factor  of 
safety*  Indeed  It  may  be  very  difficult  to  ascertain  whether 
or  not  the  readings  are  reliable,  so  the  installation  should 
be  carried  out  at  all  times  by  experienced  and  reliable 
personnel,  with  a  high  degree  of  quality  control* 


7*2  Design  sil  Instrumentation  Lead  Irsji£.b£g  for  Mis  ft  £&& 
During  the  early  stages  of  design  and  construction, 
attention  was  paid  to  this  problem  of  seepage  along  lead 
trenches*  As  a  precaution  against  this  the  addition  of 
bentonite  to  the  backfill  was  specified  and  cut-offs  were 
required  every  15  m  or  between  adjacent  piezometers, 
whichever  was  less*  The  details  of  the  actual  trench 
installations  are  given  in  Section  2*3*  These  precautions 
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may  not  have  been  adequate,  as  evidenced  by  the  behaviour  of 
piezometers  PE38*  PP36*  and  PE44*  discussed  previously  in 
Section  6*6*  Although  it  Is  unusual  practise  to  bring  lead 
trenches  into  the  upstream  shell*  perhaps  in  this  case  it 
was  rather  fortunate  as  it  provides  conclusive  evidence  of 
the  unsatisfactory  trench  behaviour*  With  this  in  mind  it 
was  decided  to  analyse  the  sensitivity  of  the  local  pore 
pressure  distributions  to  varying  trench-fill  permeability 
ra  t ios • 

Ij.2.  ggpp.ral  Hg£e-l  los.  Elgzangtgr  IcfcDcJfaeg  c 

Vaughan  (1965)  provided  an  approximate  solution  for  the 
problem  of  leakage  through  backfill  around  a  vertical 
standpipe  piezometer*  Unfortunately  this  analysis  could  not 
conveniently  be  applied  to  this  case  of  a  horizontal  pipe  in 
a  specific  flow  domain*  It  was  decided  to  model  the 
situation  using  finite  element  techniques* 

The  FPM5  computer  program  was  used  for  this  analysis* 
The  model  assumed  must  be  ax i symme t ri c *  which  presents 
problems  in  reproducing  the  asymmetric  boundary  conditions 
due  to  non— horizon ta l  flow*  The  problem  is  really  three 
dimensional  but  because  the  core  is  thin  an  assumption  of 
horizontal  flow  is  not  unreasonable*  The  cross-sectional 
area  of  the  rectangular  trench  was  approximated  by  an 
equivalent  area  circle*  The  general  model  is  shown  in  Figure 
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FIGURE  71  Instrument-Lead  Trench  Model 
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Three  basic  cases  were  studied: 

1*  The  trench  connected  to  the  upstream  shell  (CASE  1  ) 

2*  The  trench  confined  in  the  central  core  (CASE  2) 

3*  The  trench  connected  to  the  downstream  shell  (CASE  3) 
Various  permeability  and  trench  depth/core  width  ratios  were 
considered  for  each  case*  The  results  are  plotted  as 
permeability  ratio  versus  potential  values  at  specific 
locations  across  the  core*  The  effect  of  trench  diameter  is 
superimposed  on  the  charts  (Figures  7*2  to  7*4)*  Examples  of 
the  local  changes  in  pore  pressure  distribution  are  shown  on 
Figures  7*5  to  7*7* 

The  result  of  this  analysis  is  to  show  that  for  each 
case  studied  there  could  be  significant  alteration  to  the 
readings  of  piezometers  located  in  certain  critical 
positions*  The  position  is  aggravated  by  increasing  the 
trench  depth  or  permeability  ratio*  Permeability  ratios  as 
low  as  100  may  cause  significant  disturbances  to  the  pore 
pressure  regime*  This  is  consistent  with  Vaughan's  (  1965) 
conclusion  that  a  grout  to  soil  permeability  ratio  of  50  is 
satisfactory  foi4  very  low  permeability  soils  (standpipe 
installation 

These  results  are  interesting  because  as  shown  by 
Mitchell  et  al*  (1965)  (for  a  silty  clay)  it  is  relatively 
easy  to  get  a  variation  of  two  or  three  orders  of  magnitude 
in  permeability  in  the  same  material  by  simply  varying  the 
compaction  moisture  content  or  compactive  effort* 

This  analysis  shows  that  it  is  possible  to  get 
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FIGURE  7.2  Design  Chart  for  Trench  Connected  to  Upstream  Shell 
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FIGURE  7.3  Design  Chart  for  Centrally  Located  Trench 
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FIGURE  7.4  Design  Chart  for  Trench  Connected  to  Downstream  Shell 
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Typical  Finite  Element  mesh  for  Trench  Model 
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significant  local  deviations  from  the  general  pore  pressure 
distribution  within  a  core,  caused  by  the  instrument  lead 
trenches*  Therefore  it  is  necessary  to  study  individual 
trenches  in  the  Mica  dam  core  more  carefully. 


2-2.4  Modfc.l  ££e£jL£l£&lJLx  AbbIIb^  1b  U1b&  Pas  Qqjls. 

In  order  to  analyse  specific  trenches  within  the  dam 
core,  an  accurate  record  of  the  installation  method  and 
locations  of  the  lead  trenches  was  required.  This  proved 
difficult  to  obtain,  because  some  of  the  records  are  now 
more  than  10  years  old.  However  some  details  were  obtained 
from  some  unpublished  'Dam  Instrumentation  Construction 
Notes'  (Mica  dam)  from  the  years  1969  to  1972  (Caseco 
Consultants  Ltd.  1969-1972).  The  installation  procedure  was 
discussed  in  detail  in  Section  2.3,  and  layouts  of  the 
trenches  and  risers  are  shown  on  Figure  2.5. 

The  observed  behaviour  of  piezometers  PE38  and  PP36 
were  analyzed  first.  The  results  indicated  that  a 
permeability  ratio  of  about  800  is  required  to  model  the 
obse rved  behaviour.  The  resul t  for  PE 38  i s  higher  but  this 
piezometer  reads  erratically  later  on  so  it  has  been 
rejected  as  unreliable.  Next  the  model  was  applied  to 
piezometer  PE44,  once  again  the  result  indicated  a 
permeability  ratio  of  600.  This  is  very  similar  to  the 
previous  result  so  an  average  value  of  700  was  taken  for 
further  analyses.  This  infers  a  trench  permeability  of 


7xl0-4  cm/sec.See  Figure  7.8. 
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The  trench  backfill  was  sand— ben toni te— cement  as 
discussed  previously  (Section  2*3)*  In  the  upper  trenches 
the  complete  trench  was  filled  with  this  material,  Including 
cut-offs*  From  the  limited  data  available  it  appears  that  in 
many  cases  the  fill  was  compacted  very  dry  of  optimum,  and 
the  sand  used  was  likely  screened  shell  material*  It  is 
quite  possible  under  these  conditions  that  the  bentonite  did 
not  work  as  anticipated*  Rollins  and  Dylla  (1970)  and 
Mlddlebrooks  et  al*  (1978)  both  indicate  that  bentonite  may 
either  lose  its  sealing  power  or  simply  be  subject  to  piping 
under  certain  conditions*  Neither  paper  deals  specifically 
with  instrumentation  trenches  but  the  mechanism  may  still  be 
possible*  In  view  of  the  discussion  above  it  is  quite 
feasible  for  the  trenches  to  have  this  high  permeability*  It 
should  be  remembered  too  that  if  the  core  material  should 
have  a  lower  permeability  than  that  assumed  then  the  trench 
permeability  need  not  be  so  high,  because  the  disturbance  is 
caused  only  by  the  ratio  of  permeabilities* 

Next  the  lower  trenches  were  analyzed,  assuming  a 
trench  permeability  of  7  x  10-4  cm/sec*  It  was  immediately 
evident  that  the  correct  pore  pressure  trend  was  resulting* 
Particularly  for  the  cases  where  the  trenches  were  truncated 
within  the  core  (  PE25  ...  PH7;  PE36  ...  PHll;  PE37  ... 

PH12)*  For  the  cases  where  the  trench  continued  into  the 
downstream  shell  ( PE35  «•«  PH10;  PE24  **•  PH2)  the 
distribution  was  an  improvement,  particularly  on  the 
upstream  side,  however  the  downstream  results  were  too  low* 
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The  trench  PE24  •••  PH2  is  a  special  case  because  the 
permeability  of  this  material  is  assumed  to  be  lower  due  to 
the  higher  compaction  moisture  content*  Therefore  the 
permeability  ratio  is  7000  in  this  example#  The  trench  depth 
is  also  smaller  because  of  different  installation  procedure 
in  this  zone#  Some  typical  equlpotential  distributions 
resulting  from  the  analysis  for  various  trenches  are  shown 
on  Figures  7#9  to  7#11# 

A  careful  study  of  the  results,  Figures  7#12  and  7#13, 
show  a  consistant  trend#  In  the  cases  where  a  riser  exists 
below  the  trench  ( PE25  •••  PH 7 ;  PE36  •••  PH  11;  PE37  ••» 
PH12)  the  predicted  readings  were  all  a  little  high#  Where  a 
riser  is  above  the  trench  ( PE24  •••  PH2;  PE35  •••  PH10  )  the 
predicted  readings  in  the  downstream  zone  are  all  much  too 
low,  but  the  upstream  zones  are  good#  So  perhaps  the  risers 
are  influencing  the  results  also#  The  higher  permeab  ility 
being  caused  by  the  construction  technique  (Section  2#3) 
noted  around  the  risers#  This  was  much  more  difficult  to 
analyze • 

Specific  inflow  values  were  specified  at  the  nodes 
where  risers  were  located#  The  results  were  improved  in  all 
cases#  The  required  permeability  ratio  of  the  annulus  around 
the  risers  was  backcalculated  from  the  flow  value,  the  area, 
and  the  known  potentials  at  each  end  of  the  riser,  this 
yielded  reasonable  ratios  of  about  300#  These  results  are 
also  shown  on  Figures  7#12  and  7#13# 
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FIGURE  7.12  Theoretical  Pore  Pressure  Distribution  Compared  to 

Actual ,  PE24 . . .  PH  2 
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Theoretical  Pore  Pressure  Distribution  Compared  to 
Actual,  PE25.  . .  PH7,  PE  35,..  PH  10,  PE  37...  PH  12. 
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7 1  5  Conclusl ops 

Whereas  none  of  the  previous  analyses  predicted  the 
observed  pore  pressure  distribution  in  the  core,  this 
analysis  assuming  the  instrument  lead  trenches  to  have  a 
higher  permeability  than  the  surrounding  core,  provides  a 
pore  pressure  distribution  similar  to  the  observed  values* 

The  differences  are  reduced  by  including  some  effect 
from  the  risers  •  It  is  virtually  impossible  to  refine  the 
analysis  further,  because  of  the  sensitivity  of  the  results 
to  actual  construction  methods  end  variations,  details  of 
which  are  not  available*  PH2  is  the  only  piezometer  which 
does  not  fit  this  model  well*  A  way  to  explain  its  high 
readings  is  by  virtue  of  some  constructional  deviation  which 
may  not  have  appeared  significant  enough  to  report,  such  as 
backfilling  downstream  of  PH2  with  a  slightly  more 
impermeable  material* 

Additional  effects  such  as  pore  pressure  increases  due 
to  reservoir  loading  are  not  simple  to  evaluate  in  the 
vicinity  of  the  trenches,  except  that  their  magnitude  will 
probably  be  reduced*  The  anomalies  noted  in  the  transient 
seepage  analysis  (the  very  rapid  response  of  the  downstream 
piezometers)  could  very  likely  be  due  to  these  trenches  with 
a  much  higher  k/N  ratio*  The  apparent  rapid  dissipation  of 
the  end  of  construction  pore  pressures  in  the  central  core 
zone  (compared  to  the  edge  zones)  may  also  be  explained  by 
the  proximity  to  the  zones  of  higher  permeability* 
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It  Is  rather  difficult  to  provide  more  than  a 
qualitative  assessment  on  the  influence  of  the  instrument 
trenches  and  risers*  Quantitative  results  are  too  sensitive 
to  local  construction  methods  and  variations*  However  it  has 
been  adequately  demonstrated  that  the  Instrumentation 
trenches  Installed  in  the  Mica  dam  core  may  to  a  large 
extent  be  responsible  for  the  unusual  piezometric  readings* 
It  is  therefore  likely  that  the  general  pore  pressure 
distribution  within  the  core  is  considerably  different  than 
that  suggested  by  the  piezometers  and,  in  fact,  may  be  very 
similar  to  the  di stribution  suggest  ed  by  the  steady  state 
analysis* 

Several  solutions  are  possible  to  reduce  the  likelihood 
of  pref ferential  seepage  paths  along  the  instrument-lead 
trenc he  s • 

-  The  backfill  material  for  the  trenches  should  be 
the  same  as  the  surrounding  core  material,  with  the 
larger  soil  sizes  removed*  The  compaction  moisture 
content  and  density  should  be  maintained  as  close  as 
possible  to  the  surrounding  core  material*  This  will 
aid  in  reducing  hydraulic  property  differences* 

Strict  construction  supervision  is  required  at  all 
times  and  deviation  from  the  specifications  should 
not  be  allowed  without  careful  consideration  of  the 
possible  consequences* 

—  In  very  narrow  cores  of  low  permeability  special 
precautions  may  be  necessary  to  ensure  Instrument 
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reliability.  Procedures  such  as  those  suggested  by 
Blight  (1970)  may  be  satisfactory. 

—  For  cases  where  instrumented  sections  are  not  far 
apart  along  the  crestt  the  instrument-lead  trenches 
could  run  parallel  to  the  crest.  Special  attention 
to  detail  would  be  required  where  the  leads  turn 
perpendicular  to  the  crest  and  exit  into  the 
downstream  shell.  This  alternative  is  controlled  by 
the  economics  of  the  situation. 
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Significance  of  the  Processes  I_iiYfiI.3tc.ji  during  Resevoi r 

Ijgfififinctwfiflt  ansi  Iran 9 lent  Safifi&ast*. 

-gj-l-r.l  Xfilal  stress  changes . 

The  total  stress  changes  along  the  upstream  face  of  the 
core  caused  by  reservoir  filling  will  be  less  than 
hydrostatic  values*  These  stress  increments  will  cause  a 
change  in  pore  pressure  distribution  throughout  the  core 
which  is  significant  but  apparently  does  not  dominate  the 
observed  pore  pressure  distribution*  The  estimated  values 
are  functions  of  the  induced  shear  stress  distribution 
assumed  across  the  core*  In  this  thesis  the  assumed 
distribution  was  compared  to  the  analysed  case  of  the 
Orville  Dam  and  the  measured  case  of  the  Gespatsch  Dam* 


gj.lj.2  seesam 

The  transient  seepage  phase  was  estimated  to  be 


relatively  rapid  in  the  Mica  dam  c a se ,  taking  approximately 
3  years  despite  the  low  permeability  of  the  core  (10  — 7 
cm/sec)*  The  advance  of  the  phreatic  surface  corresponds 
reasonably  closely  to  the  observed  water  levels  in  the 
movement  gauges*  The  estimated  steady  state  phreatic  surface 
by  conventional  analysis  is  very  much  higher  in  the 
downstream  zone  of  the  core  when  compared  to  an  analysis 
which  considers  flow  within  the  capilliary  fringe  (Freeze 
1971).  This  anomaly  is  most  pronounced  in  a  dam 
configuration  similar  to  the  Mica  dam*  The  downstream 


movement  gauge  appears  to  substantiate  the  argument 


for  a 
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lower  phreatic  surface  and  therefore  the  transient  analysis 
must  be  considered  inaccurate  near  the  downstream  core  face* 


8*1*3 


JErftctmringt 


It  has  been  estimated  that  alt 
fracturing  is  feasible  upon  reservo 
in  total  stresses  are  of  sufficient 
phenomenon  unlikely*  Hydraulic  free 
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XQ-H  glfflBli.tj-ag 


For  the  purposes  of  this  study  the  processes  leading  to 
pore  pressure  changes  during  reservoir  filling  and  transient 
flow  were  separated  and  analyzed  individually  with 
simplifications*  There  are  many  ill— defined  parameters  such 
as  the  storage  coefficient,  path  dependent  pore  pressure 
parameters,  which  preclude  the  use  of  more  refined  finite 
element  analysis*  It  would  be  of  limited  benefit  to  refine 
the  analysis  until  the  significant  mechanisms  and  parameters 
are  better  defined  and  understood*  It  may  be  very  difficult 
to  model  the  behaviour  of  a  core  in  the  laboratory  in  order 
to  better  determine  the  relevant  parameters  because  the 
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results  are  very  sensitive  to  in— situ  conditions  such  as 
conpactlon  effort,  compaction  moisture  content, 
permeability,  porosity  and  compressabl lity ,  some  of  which 
change  with  time  depending  on  construction  history  and 
location  in  the  embankment* 

Suitable,  well  placed  instrumentation  would  help  to 
provide  further  insight  into  the  contributions  of  each 
seperate  mechanism* 


£U 2  SagRftrlgflfl  gilh  Observed  Behaviour 

It  is  apparent  that  the  observed  pore  pressure 
distribution  within  the  Mica  dam  core  is  at  variance  with 
the  calculated  steady  state  distribution,  even  though  the 
evidence  suggests  that  essentially  steady  state  conditions 
have  been  achieved  after  about  3  years* 

The  observed  distribution  is  not  reflected  by  the 
reservoir  loading  nor  the  transient  seepage  analysis*  It 
must  be  concluded  that  some  alternative  phenomenon  was 
producing  the  observed  distribution*  A  study  was  completed 
which  indicates  the  sensitivity  of  local  pore  pressures  to 
non— homogeneities  present  in  the  form  of  i n s tr ument— lead 
trenches*  The  evidence  from  the  piezometers  connected  to  the 
upstream  shell  shows  that  these  trenches  must  provide  a 
preferential  seepage  path*  This  model  when  applied  to 
specific  trenches  within  the  till  core  readily  reflected  the 
observed  trend*  The  magnitudes  are  sensitive  to  permeability 
ratios  and  specific  trench  sizes,  accurate  details  of  which 
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are  very  difficult  to  obtain* 

Similarly  if  the  risers  are  considered  to  behave  as 
preferential  seepage  paths  the  agreement  between  the 
calculated  steady  state  pore  pressure  distribution  and  the 
measured  distribution  is  improved*  It  would  be  useless  to 
refine  the  analysis  further  because  the  specfic  construction 
details  required  for  input  are  not  available*  The  analysis 
indicates  in  a  quantitative  manner  that  the  piezometers 
placed  in  the  Mica  dam  core  are  not  reflecting  the  general 
pore  pressure  distribution  during  impounding*  The 
non— homogeneities  caused  by  the  instrument-lead  trenches  are 
sufficient  to  create  a  local  pore  pressure  disturbance, 
which  is  being  monitored  by  the  piezometers*  It  is  concluded 
that  the  average  pore  pressure  distribution  within  the  core 
is  probably  similar  to  that  proposed  by  steady  state  theory, 
with  possible  errors  in  the  downstream  core  zone  depending 
on  whether  the  flow  in  the  capilliary  fringe  is  being 
consi dered • 

This  emphasizes  the  very  considerable  care  required 
when  installing  piezometers  in  a  narrow  core  embankment*  If 
sufficient  protection  is  not  employed  against  preferential 
seepage  the  results  will  likely  be  valueless  and  the  expense 
of  the  entire  installation  wasted* 

A  serious  consequence  of  inadequate  installation 
procedure  would  be  the  increased  likelihood  of  hydraulic 
fracture  Initiating  if  the  trench  were  connected  to  the 
upstream  shell*  Therefore  designers  should  never  consider 
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running  "trenches  from  "the  lower  core  info  "the  upstream 
shell,  and  this  should  not  be  done  in  the  upper  core  without 
considerable  Justification  and  consideration  of  the  possible 
consequences#  Considerable  attention  is  required  during  the 
design  and  installation  to  ensure  that  the  hydraulic 
properties  of  the  1 nst rumen t— lead  trenches  are  similar  to 
the  surrounding  core  material* 


8-f  3  Sgl.ecllgQ  OJL  Piezometer  Type 

Some  conclusions  are  drawn  regarding  piezometer  types 
used  in  embankment  dams* 

Hydraulic  *  embankment  type*  piezometers  are  excellent 
and  should  be  used  where  possible*  Their  main  attraction  is 
simplicity  of  operation,  the  facility  for  de- airing  and 
in— situ  testing*  If  adequate  care  is  taken  during 
installation  their  performance  Is  completely  satisfactory. 

Electric  piezometers  require  special  attention  in 
partly  saturated  zones*  The  choice  of  a  high  air  entry  stone 
does  not  remove  the  possibility  that  the  tip  is  recording 
pore  air  pressure  in  zones  of  low  stress*  Their  design 
precludes  the  possibility  of  flushing  or  in— situ  testing* 
They  have  been  proven  susceptible  to  overvoltage  failures 
from  lightning  strikes  unless  very  special  attention  is  paid 
to  overvoltage  protection* 

Pneumatic  piezometers  generally  lack  a  deairing 
facility  and  in— situ  tests  are  normally  impossible*  Long 
lead  lengths  and  certain  lead  types  may  present  problems. 
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Operator  error  is  not  uncommon*  Insufficient  evidence  is 
available  in  the  literature  regarding  the  long  term 
behaviour  and  accuracy  of  these  piezometers* 

In  many  instances  electric  and  pneumatic  piezometers 
may  be  excellent  choices  but  it  is  the  author's  opinion  that 
where  practical  hydraulic  'embankment  type'  piezometers 
should  be  used  in  partly  saturated  compacted-fill 
embankments*  Their  suitability  and  long  term  reliability 
have  been  adequately  demonstrated  in  the  literature*  If  a 
variety  of  piezometer  types  are  employed  then  care  must  be 
exercised  to  ensure  that  different  piezometer  types  are 
placed  within  zones  of  similar  potential*  In  a  few  locations 
it  is  desirable  to  provide  a  hydraulic  piezometer  close  to 
the  alternate  type,  to  serve  as  a  calibration  check*  In  this 
manner  an  adequate  description  of  the  pore  pressure 
distribution  may  still  be  possible  If  a  particular 
piezometer  type  malfunctions*  It  must  be  remembered  that 
once  a  malfunction  occurs  it  may  be  impossible  to  replace 
that  piezometer  and  it  will  certainly  be  very  expensive  to 
do  so • 

8.4  Recommendations  for  F\»rth£r  &££L£&XL£l2 

Based  on  the  work  presented  here  the  following  is 
recommended  for  further  research* 

1*  The  effect  of  rotation  of  principal  stresses  on 
the  pore  pressure  parameters  require  further 


attention • 
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2*  Further  Investigation  is  required  into  the  shape 
of  the  steady  state  phreatic  surface  in  a  core 
where  the  unsaturated  or  capillary  fringe  flow 
is  considered*  Particularly  useful  would  be 
investigation  of  case  histories  which  indicate 
the  phreatic  surface  with  time* 

3*  In  order  that  transient  flow  may  be  succesfully 
analysed  it  is  necessary  to  provide  a  more 
accurate  description  of  the  storage  coefficient 
and  its  sensitivity  to  parameters  such  as  degree 
of  saturation,  compation  water  content, 
porosity,  degree  of  compaction,  stress  history, 
and  grain  size  distribution* 

4*  A  particularly  important  area  of  research  would 
be  an  Investigation  of  hydraulic  fracture.  It 
may  be  possible  to  use  induced  hydraulic 
fracture  as  an  indicator  of  degree  of  arching  in 
critical  zones  of  the  core*  This  in  turn  could 
be  related  to  the  likelihood  of  hydraulic 
fracture  occurring  on  reservoir  filling  •  Earth 
pressure  cells  have  not  been  very  sucessful  in 
embankments  to  date,  but  perhaps  they  could  be 


u 


sed  in  conjunction  with  hydraulic  fracture 


tests  to  determine  in— situ  stress  states  more 


accurately 


- 
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